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ABSTRACT 

The zodiacal cloud is a thick circumsolar disk of small debris particles pro- 
duced by asteroid collisions and comets. Their relative contribution and how 
particles of different sizes dynamically evolve to produce the observed phenom- 
ena of light scattering, thermal emission, and meteoroid impacts are unknown. 
Until now, zodiacal cloud models have been phenomenological in nature, com- 
posed of ad-hoc components with properties not understood from basic physical 
processes. Here, we present a zodiacal cloud model based on the orbital proper- 
ties and lifetimes of comets and asteroids, and on the dynamical evolution of dust 
after ejection. The model is quantitatively constrained by IRAS observations of 
thermal emission, but also qualitatively consistent with other zodiacal cloud ob- 
servations, with meteor observations, with spacecraft impact experiments, and 
with properties of recovered micrometeorites. We find that 85-95% of the ob- 
served mid-infrared emission is produced by particles from the Jupiter-family 
comets (JFCs) and <10% by dust from long period comets. The JFC parti- 
cles that contribute to the observed cross-section area of the zodiacal cloud are 
typically D ^ 100 /xm in diameter. Asteroidal dust is found to be present at 
<10%. We suggest that spontaneous disruptions of JFCs, rather than the usual 
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cometary activity driven by sublimating volatiles, is the main mechanism that 
hbrates cometary particles into the zodiacal cloud. The ejected mm to cm-sized 
particles, which may constitute the basic grain size in comets, are disrupted on 
<10,000 yrs to produce the 10-1000 fim grains that dominate the thermal emis- 
sion and mass influx. Breakup products with D > 100 /im undergo a further 
collisional cascade with smaller fragments being progressively more affected by 
Poynting-Robertson (PR) drag. Upon reaching D < 100 /xm, the particles typ- 
ically drift down to <1 AU without suffering further disruptions. The resulting 
Earth impact speed and direction of JFC particles is a strong function of par- 
ticle size. While 300 //m to 1 mm sporadic meteoroids are still on eccentric 
JFC-like orbits and impact from antihelion/helion directions, which is consistent 
with the aperture radar observations, the 10-300 fim particles have their orbits 
circularized by PR drag, impact at low speeds and are not detected by radar. 
Our results imply that JFC particles represent ~85% of the total mass influx 
at Earth. Since their atmospheric entry speeds are typically low (^14.5 km s"-^ 
mean for D — 100-200 //m with f=il2 km s~^ being the most common case), many 
JFC grains should survive frictional heating and land on the Earth's surface. 
This explains why most micrometeorites collected in antarctic ice have primi- 
tive carbonaccoTis composition. The present mass of the inner zodiacal cloud 
at <5 AU is estimated to be 1-2 x 10^^ g, mainly in D = 100-200 /xm parti- 
cles. The inner zodiacal cloud should have been >10^ times brighter during the 
Late Heavy Bombardment (LHB) epoch ^3.8 Gyr ago, when the outer planets 
scattered numerous comets into the inner solar system. The bright debris disks 
with a large 24-/xm excess observed around mature stars may be an indication 
of massive cometary populations existing in those systems. We estimate that 
~10^^, ~2 X 10^^ and ~2 x 10^° g of primitive dark dust material could have 
been accreted during LHB by the Earth, Mars and Moon, respectively. 

Subject headings: Zodiacal light; Comets: dust; Debris disks; Meteorites 



1. Introduction 

The zodiacal cloud is a dynamic assembly of meteoroids in bound orbits around the 
Sun. The orbits depend on particle size, location in the cloud, and the type of parent body. 
Interstellar dust particles that pass through the solar system are not considered in this paper, 
nor are small meteoroid fragments that move out of the solar system on hyperbolic orbits 
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( "beta-meteoroids" ) . 

Traditionally, the zodiacal cloud has been described with phenomenological models of 
dust distributions to explain the amount of scattered light (Hong, 1985; Kniessel and Mann, 
1991; Ishiguro et al. 1999; Hahn et al., 2002), the Doppler shifts of the solar Mg I Fraunhofer 
line (Hirschi and Beard, 1987, Mukai and Mann, 1993; Clarke et al., 1996; Reynolds et al., 
2004), and the more easily to interpret thermal emission observed in various lines of sight 
(Kelsall et al., 1998; Maris et al., 2006). Particularly good scattered light observations came 
from the Clementine mission (Hahn et al., 2002), while thermal infrared observations are 
mostly from the Infrared Astronomical Satellite - IRAS (Low et al., 1984; Hauser et al., 1984; 
Good et al., 1986; Sykes et al., 1990), the Cosmic Background Explorer - COBE (Reach et 
al., 1995; Kelsall et al., 1998), the Midcourse Space Experiment - MSX (Price et al., 2003), 
the Infrared Space Observatory - ISO (Fixsen and Dwek, 2002; Leinert et al., 2002; Reach et 
al., 2003; Mueller et al. 2005), and the Spitzer Space Telescope (Bhattacharya and Reach, 
2004; Reach et al., 2007). 

The femenological models successfully describe the size, spatial and velocity distribu- 
tions of dust particles in the solar system (e.g., Griin et al., 1985; Divine, 1993; Dikarev et 
al., 2004). They are particularly useful for accessing the satellite impact hazard, designing 
spacecraft impact experiments and studies of extrasolar emission sources such as the cosmic 
microwave background (e.g., Kelsall et al., 1998). The femenological models, however, fall 
short in answering basic questions related to the origin of the zodiacal cloud, its temporal 
brightness variability, and the provenance of interplanetary particles collected at the Earth. 
Consequently, the origin of the zodiacal cloud, interplanetary dust particles (IDPs) collected 
in the Earth's stratosphere (e.g.. Love and Brownlee, 1994) and micrometeorites (MMs) on 
the ground (Taylor et al., 1996; Engrand and Maurette, 1998; Farley et al., 1998, 2006; 
Genge, 2006) is still being a matter of considerable debate. This limits our ability to link the 
detailed laboratory studies of IDPs and MMs to the properties of their parent bodies, and 
to use the zodiacal cloud as a valuable reference for studies of the exozodiacal debris disks. 

Detailed dynamical models can be more useful in this context. At the root of dynamical 
models are the physical properties of interplanetary dust, such as density, geometric albedo, 
elemental composition, mineralogy, tensile strengh, heat capacity, etc. (e.g., Dumont and 
Levasseur-Regourd, 1988; McDonnell and Gardner, 1998; Gustafson, 1994; Gustafson et 
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al., 2001; Levasseur-Regourd et al., 2001), which determine the behavior of particles in 
interplanetary space (e.g., planetary perturbations, collisions, sublimation, sputtering) and 
their interaction with a detector (e.g., ablation of micrometeorites in the Earth's atmosphere, 
^He retention, thermal radiation, light scattering). In dynamical models, the individual 
particles are tracked by numerical codes as they evolve by various processes from their 
sources (assumed to be, e.g., asteroids, comets, satellites or Kuiper belt objects) to sinks 
(e.g., when they sublimate, disrupt, impact or leave the solar system). Insights into the 
origin of the zodiacal cloud can be obtained by calibrating the results of dynamical models 
on observations. 

Until now, detailed dynamical models have been only developed for asteroidal dust to 
explain the origin of the zodiacal dust bands (e.g., Dermott et al. 1984; Grogan et al., 1997, 
2001; Reach et al, 1997; Nesvorny et al. 2006; Vokrouhlicky et al., 2008) and trapped dust 
in Earth's Langrange points first seen in IRAS observations (Dermott et al., 1994a). It has 
been established that the dust bands originate from the youngest asteroid famihes (Nesvorny 
et al., 2003, 2008). However, claims that asteroids are a major if not dominant source of 
zodiacal dust, by assuming that all main belt asteroids contribute dust (e.g., Dermott et al. 
1995; Durda and Dermott, 1997; Kortenkamp and Dermott, 1998) has remained in doubt. 

Models of the zodiacal cloud need not only explain line-of-sight properties, but also the 
observed influx of meteors (see Ceplecha et al., 1998; Jenniskens, 2006, for a review) and 
the impact rate of meteoroids on sateUites (Love and Brownlee, 1993). Until now, models 
that were developed to explain these dynamical phenomena (e.g., Griin et al., 1985; Divine 
et al, 1993; Kessler et al. 1994; Staubach et al, 1997; Dikarev et al, 2004) were based 
on ad-hoc populations of meteoroids in various types of orbits without a dynamical under- 
pinning to sources and sinks. Moreover, all current satellite impact models use meteoroid 
velocity distributions (both magnitude and spatial direction) derived from meteor observa- 
tions (Taylor and Elford, 1998; Jones and Brown, 1993; Brown and Jones, 1999; Brown and 
Campbell-Brown, 2003), which pertain to much bigger particles than typically encountered 
by sateUites. 

In this paper, we investigate what fraction of the zodiacal cloud is due to cometary ver- 
sus asteroidal dust by calculating the evolution of dust particles under solar radiation forces 
and planetary perturbations (including resonances and close encounters), ejected from model 
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populations of all potential sources (not just representative examples). The source popula- 
tions include asteroids, active and mostly dormant Jupiter- family comets (JFCs), Halley-type 
comets and long-period Oort-cloud comets. In recent years, much insight was gained into the 
dynamical characteristics of these populations and the number of asteroids and comets that 
can contribute dust to the zodiacal cloud (e.g., Levison and Duncan, 1994, 1997; Jedicke and 
Metcalfe, 1998; Wiegert and Tremaine, 1999; Dones et al., 2004; Francis, 2005; Gladman et 
al. 2009). At the same time, it was realised that mostly dormant JFCs are the main source of 
meteoroid streams in the inner solar system (Jcnniskens 2006; 2008) and responsible for the 
antihelion/helion sources in the sporadic meteoroid background (Jenniskens, 2006; Wiegert 
et al., 2009). 

Here we couple these new insights to dynamical behaviour of cometary and asteroidal 
dust in order to evaluate the contribution of various sources to the thermal emission of 
zodiacal dust and the influx of micrometeorites. We show that JFCs are the main source 
of zodiacal dust inside 5 AU and the most likely source of the micrometeorites found on 
Earth. Our results also provide quantitative constraints on dust lifetimes, influx rates, and 
velocity distributions directly from the known abundances of meteoroid parent bodies. The 
results are used to quantify the properties of zodiacal dust cloud in the past and discuss 
imphcations for studies of exozodiacal debris disks. 

To set up the stage for our modeling described in §3, we discuss IRAS observations 
of the zodiacal cloud in §2. Results are presented in §4. In §5, we estimate the current 
and historical terrestrial accretion rates of dust and discuss the implications of our work for 
studies of micrometeorites and debris disks. Comet disruptions/splitting events are reviewed 
in §6. We suggest that they are the main mechanism by which the cometary particles are 
liberated from their parent bodies into the zodiacal cloud. Previous work and origin of dust 
particles beyond Jupiter are discussed in §7 and §8, respectively. 

2. Constraints 

Our primary constraints arc the observations of the zodiacal cloud by IRAS which have 
been confirmed by COBE and Spitzer (Hauscr ct al., 1984; Low et al., 1984; Kelsall et al., 
1998; Sykes et al., 2004). IRAS measured mid-infrared (MIR) fluxes in four filters with 
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effective wavelengths of 12, 25, 60 and 100 /im. These filters can be used as windows into 
the dust distribution at different distances from the Sun. Measurements in the 12-/im IRAS 
band are mainly sensitive to distributions of particles at ~l-2 AU, while the 25- and 60- 
Hm band measurements preferentially detect thermal emission from larger distances. IRAS 
observations in the 100- /xm band are less useful for probing the thermal radiation of dust 
particles in the inner solar system. 

IRAS showed that the MIR brightness of the zodiacal cloud peaks at the ecliptic and 
has broad wings that extend all the way to the ecliptic poles (Fig. [T]). The variation with the 
ecliptic longitude is minimal indicating that the zodiacal cloud is a nearly symmetrical disk 
of circumsolar particles that is roughly centered at the Sun (e.g., Staubach et al., 2001). The 
significant flux received from the ecliptic poles (~l/3-l/4 of the ecliptic flux) also suggests 
that the cloud must be rather thick in the normal direction to the ecliptic plane. 

Following Nesvorny et al. (2006; hereafter NVBS06), we selected several representative 
IRAS scans that met the following conditions. (1) We used scans that covered a continuous 
range of ecliptic latitudes from b < —80° to 6 > 80°. (2) We did not use scans that had 
gaps created when the telescope skipped over bright sources. (3) We did not use scans that 
showed strong emission from extrasolar sources such as the galactic plane, galactic cirrus, 
point sources, etc. (4) We required that the selected scans covered all values of ecliptic 
longitude and the available range of solar elongations, /©. Tables 2 and 3 in NVBS06 list 
the basic information about the selected scans. 

In this work, we only consider scans with Iq ^ 90°. This is mainly done because the 
detail variation of brightness with Iq is difficult to characterize as it requires an appropriate 
model for the collisional disruption of particles (NVBS06). We do not include the effects of 
collisions between particles in our model (except in §4.2, where a simple model for collisions 
is used). The considered value of Iq is in the middle of the available range and thus best 
represents IRAS observations. 

The zodiacal cloud is known to be warped and have a center that is slightly offset from 
the Sun. These features are produced by gravitational perturbations from Jupiter (Dermott 
et al., 1995, 2001). Since we are not interested in these detailed features of the zodiacal cloud 
here, we removed them by combining the selected IRAS scans into a representative profile. 
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This was done by first shifting tlie scans by a small value in latitude (<2°), so that the peak 
of emission was centered exactly at 6 = 0, and calculating the average flux from all selected 
scans as a function of b. We ended up with profiles showing the mean fluxes at 12, 25 and 
60 fim wavelengths as a function of ecliptic latitude (Fig. [2]). These profiles represent the 
main constraints on the work described here. Additional constraints are discussed in §4.3. 

3. Model 

Our model of the zodiacal cloud has four parts: we (i) define the initial orbital distribu- 
tions of particles for different sources (asteroid and comet populations); (ii) track the orbital 
evolution of particles with various sizes from sources to sinks; (iii) determine the thermal 
infrared emission from these synthetic particle distributions; and (iv) model the detection of 
their emission by IRAS. These model components are described below. 

To simplify things, we do not initially account for the collisional disruptions of dust parti- 
cles in the model. This is a major assumption which we verify in §4.2. For example, NVBS06 
showed how the collisional disruption of particles, and production of smaller daughter prod- 
ucts, affect the spatial distribution of particle populations in the asteroidal dust bands. The 
work described here, however, is less sensitive to collisional processes because the overall 
gross shape of the zodiacal cloud should be mainly controlled by the orbital distribution of 
its source population(s), violent dynamics of particles moving on planet-crossing orbits, and 
the effects of PR drag. These features/processes are included in our model as we describe 
below. 

3.1. Sources 

Our model starts with the source population of objects. This may be either of the 
following: (1) individual asteroid groups such as the Karin, Veritas or Beagle asteroid families 
(NVBS06, Nesvorny et al., 2008); (2) asteroid belt as a whole; (3) active JFCs defined by 
their debiased orbital distribution obtained from Levison & Duncan (1997; hereafter LD97) 
and their physical lifetime; (4) JFCs with orbital distribution similar to (3) but dynamically 
evolved beyond their nominal active lifespan; (5) Halley-type comets (HTCs); and (6) Oort- 
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cloud comets (OCCs). These source populations are described in a more detail below. 

We ignore the contribution of interstellar dust particles because thermal emission from 
these small particles (diameter D < 1 fim) would create diagnostic spectral features in the 
mid-infrared wavelengths that are not observed (e.g., Reach et al., 2003). More specifically, 
observations of the MIR spectrum of the zodiacal cloud by Reach et al. (2003) suggest that 
the bulk of the zodiacal cloud is produced by ~1 0-100 /im particles. We also ignore dust 
produced by disruptive collisions in the Kuiper belt (hereafter KB dust) because Landgraf 
et al. (2002) and Moro-Martin and Malhotra (2003) showed that KB dust particles should 
represent only a minor contribution to the inner zodiacal cloud. 

Results for (1) were taken from NVBS06 and Nesvorny et al. (2008) who found that 
the three main dust bands discovered by Low et al. (1984) originate from the Karin, Veritas 
and Beagle asteroid families. According to these results (Fig. [1]), the three main dust bands 
represent only ^9-15% of the zodiacal cloud emission at the ecliptic and <5% overall. About 
a dozen other dust bands have been identified (Sykes, 1990). Since these dust bands are much 
fainter than the three main dust bands, Nesvorny et al. results set the upper limit on the 
contribution of identified asteroid breakups to the zodiacal cloud. 

The orbital distribution of the asteroid belt source (2) is modeled by using the observed 
orbital distribution of asteroids with D > 15 km. We obtained this distribution from the 
ASTORB catalog (Bowell et al., 1994). This sample should be complete and unbiased 
(Jedicke and Metcalfe, 1998; Gladman et al., 2009). The orbital distribution of asteroids 
with D < 15 km, which may be a better proxy for the initial distribution of dust produced in 
main belt collisions (Sykes and Greenberg, 1996), is roughly similar to that of large asteroids. 
Thus, the orbital distribution that we use should be a reasonable assumption for (2). 

The orbital distributions obtained by numerical integrations of test particle trajectories 
from individual comets would suffer from the heavily biased comet catalogs. Moreover, there 
are hundreds of known comets, each producing dust at a variable and typically unknown rate. 
In addition, it is also possible that the present zodiacal dust complex contains particles from 
lost parents as suggested by the orphaned Type-II trails (Sykes, 1990) and identification 
of meteoroid streams with disrupted JFCs (see Jenniskens (2008) for a review). In view of 
these difficulties, we resorted to the following strategy. 



- 9 - 



The orbital distribution of JFCs was taken from LD97 who followed the evolution of 
bodies originating in the Kuiper belt as they are scattered by planets and evolve in small 
fractions into the inner solar system. Starting at the time when the comets' perihelion 
distance, g, first drops below 2.5 AU in the LD97 simulations, we include it as an active 
JFC in our list of source objects. LD97 showed that the orbital distribution of visible 
JFCs obtained in this way nicely approximates the observed distribution. Moreover, LD97 
argued that the inclination distribution of new JFCs (reaching q < 2.5 AU for the first time) 
is relatively narrow. The inclination distribution widens at later times as the JFC orbits 
become more spread by Jupiter encounters. 

By comparing the width of the model inclination distribution with that of the observed 
JFCs, LD97 were able to estimate the fading lifetime of JFCs, tjpc; defined as the char- 
acteristic time between their first and last apparitions. They found that tjFc = 12, 000 yr 
with a 90% confidence interval 3, 000 < tjpc < 30, 000 yr. See Fig. [3] for an illustration of 
the steady-state JFC orbit distribution for tjFc = 12, 000 yr. This is our initial distribution 
of particles produced by active JFCs. We use tjFC a free parameter in our model with 
values extending to tjpc = 100, 000 yr to account for the possibility that an important dust 
component may be produced by old dormant JFCs as they spontaneously disrupt. 

Our models for the orbital distributions of HTCs and OCCs are simpler than the one 
described for JFCs above, because we do not have in hands an appropriate numerical model 
that we could use with confidence. Fortunately, this is not a major limitation factor in this 
work because our main results described in §4 are not sensitive to the detailed properties of 
the HTCs and OCCs populations. 

For HTCs, we assume that the differential distribution of the perihelion distance, dA^(g), 
is dN{q) oc qdq, and set an upper limit of g at 3 AU. HTCs typically become visible/active 
only if they reach g < 3 AU. Similarly, the cumulative semimajor axis distribution of HTCs, 
A^(< a), is assumed to be A^(< a) oc a with an upper cut on a at 50 AU. The differential 
inclination distribution dN{i) is taken from Levison et al. (2006). This distribution can be 
approximately described by 

dN{i) oc sin(i)e~^(^) di (1) 

with a = 30°. 
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According to Francis (2005), the long-period comets have dN{q) oc (1 + y^) d^ for 
q < 2 AU. For q > 2 AU, Francis' study predicts dN{q) being flat or dechning while we 
would expect the perihelion distribution to increase with q. It probably just shows that the 
distribution is not well constrained for > 2 AU. We use dN{q) oc 2Al{q/2)'^ dq ior q > 2 AU 
with < 7 < 1. The semimajor axis values of OCCs are set between 10,000 and 50,000 AU, 
which is known as the Oort spike (Wiegert and Tremaine, 1999). We also use a— 1, 000 AU 
to check on the dynamical behavior of dust particles launched from the returning OCCs. 
The inclination distribution of OCCs is set to be uniformly random between and 180°. 

We use two different methods to launch particles from their source objects. In the 
flrst method, chosen to approximate the ejection of particles from active comets, we launch 
particles at the perihelion when the mean anomaly M — 0. In the second method, we launch 
particles along the orbit with uniform distribution in M. This second method is more 
appropriate for the asteroidal particles and for particles produced by comet disruptions. 
Indeed, the identifled comet disruptions do not seem to be correlated in any way with 
the perihehon passage (e.g., Weissman, 1980). To simphfy things, we neglect the ejection 
velocities of dust particles from their parent object and assume that they will initially follow 
the parent object's orbit modifled by the radiation pressure. 

The individual comets in our model are assumed to contribute in roughly the same 
proportion to the circumsolar dust complex. The reasoning behind this assumption is that if 
an individual super-comet were the dominant source of circumsolar dust, the zodiacal cloud 
would not have such a smooth and symmetrical structure. In fact, the observed smooth 
structure of the zodiacal cloud probably implies a source population that contains numerous 
objects that are well mixed in the orbital space. 

3.2. Orbit Evolution 

The orbits of particles were tracked using the Wisdom-Holman map (Wisdom and Hol- 
man, 1991) modified to include effects of radiation forces (Burns et al., 1979; Bertotti et al., 
2003). The acceleration F on a particle due to these forces is 




V 



c 



(2) 
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where R is the orbital radius vector of the particle, V is its velocity, G is the gravitational 
constant, is the mass of the Sun, c is the speed of light, and R = dR/ dt. The acceleration 
([2]) consists of the radiation pressure and velocity-dependent PR terms. Parameter (3 is 
related to the radiation pressure coefficient, Qpr, by 

= 5.7 X 10-^ ^ , (3) 
ps 

where radius s and density p of the particle are in cgs units. Pressure coefficient Qpj- can be 
determined using the Mie theory (Burns et al., 1979). We used Qpr = 1 which corresponds 
to the geometrical optics limit where s is much larger than the incident-light wavelength. 
We assumed that the solar-wind drag force has the same functional form as the PR term 
and contributes by 30% to the total drag intensity (Gustafson, 1994). 

We used particles with diameter D = 2s = 10, 30, 100, 200, 300, 1000 pm and set 
their bulk density to p = 2 g cm~^. For comparison. Love et al. (1994) reported p ^ 2 
g cm~^ for stratospheric-collected IDPs, while McDonnel and Gardner (1998) found mean 
p = 2-2.4 g cm~^ from the analysis of data collected by the LDEF and Eureca satellites. 
On the other hand, density of 1 g cm~^ has been often assumed for cometary matter (e.g., 
Joswiak et al, 2007; Wiegert et al., 2009). Griin et al. (1985) suggested that 20-40% of 
particles may have low densities whereas most meteoroids have p = 2-3 g cm^^. Our results 
may be easily rescaled to any p value and we explicitly discuss the effect of p wherever it is 
appropriate. 

The particle orbits were numerically integrated with the swif t_rmvs3 code (Levison and 
Duncan, 1994) which is an efficient implementation of the Wisdom-Holman map and which, 
in addition, can deal with close encounters between particles and planets. The radiation 
pressure and drag forces were inserted into the Keplerian and kick parts of the integrator, 
respectively. The change to the Keplerian part was trivially done by substituting by 
niQ^l — P), where /3 is given by Eq. ([3]). 

The code tracks the orbital evolution of a particle that revolves around the Sun and 
is subject to the gravitational perturbations of seven planets (Venus to Neptune) until the 
particle impacts a planet, is ejected from the solar system or drifts to within 0.03 AU from 
the Sun. We removed particles that evolved to R < 0.03 AU because the orbital period for 
R < 0.03 AU is not properly resolved by our integration timestep. In 4.2, we also consider 



the effect of coUisional disruption by removing particles from the simulations when they 
reach their assumed physical hfespan. We followed 1,000-5,000 particles for each D, source, 
and parameter value(s) that define that source. 



3.3. Thermal Emission of Particles 

Particles were assumed to be isothermal, rapidly rotating spheres. The absorption was 
assumed to occur into an effective cross-section tts^, and emission out of Atts^. The infrared 
flux density (per wavelength interval dA) per unit surface area at distance r from a thermally 
radiating particle with radius s is 

F{X) = e{X,s)B{X,Tf^, (4) 

where e is the emissivity and B{X, T) is the energy flux at (A, A -|- dA) per surface area from 
a black body at temperature T: 

B{\T) = ^-^[e^^'>^'^-l]-' . (5) 

In this equation, h = 6.6262 x lO"^'' J s is the Planck constant, c = 2.99792458 x 10^ m s"^ is 
the speed of hght, and k — 1.3807 x 10~^^ J is the Boltzmann constant. We used e—1 
which should be roughly appropriate for the large particles used in this work. See NVBS06 
for a more precise treatment of e(A, s) for dust grains composed of different materials. 

To determine the temperature of a particle at distance R from the Sun, wc used the 
temperature variations with R that were proposed by different authors from spectral observa- 
tions of the zodiacal cloud (e.g., Dumont ct al., 1988; Rcnard et al., 1995; Leinert et al., 2002; 
Reach et al., 2003). For example, Leinert ct al. (2002) proposed that T{R) = 280 /R^-^^ K 
near = 1 AU from ISOPHOT spectra. We used T{R) = T^av/R^ K, where Tiau is the 
temperature at 1 AU and 5 is a power index. We varied Tiau and S to see how our results 
depend on these parameters. Values of Tiau ~ 280 K and S — 0.5 correspond to the equi- 
librium temperature of large dark particles. Values 5 < 0.5 would be expected, for example, 
for fluffy particles with small packing factors (e.g., Gustafson et al., 2001). 
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3.4. Synthetic Observations 



To compare our results with IRAS observations described in §2, we developed a code 
that models thermal emission from distributions of orbitally evolving particles and produces 
infrared fluxes that a space-borne telescope would detect depending on its location, pointing 
direction and wavelength. See NVBS06 for a detailed description of the code. 

In brief, we define the brightness integral along the line of sight of an infrared telescope 
(defined by fixed longitude / and latitude b of the pointing direction) as: 



where r is the distance from the telescope, F{D,r) is the infrared flux (integrated over 
the wavelength range of the telescope's system) per unit surface area at distance r from 
a thermally radiating particle with diameter D. S{R,L,B) defines the spatial density of 
particles as a function of the heliocentric distance, R, longitude, L, and latitude, B (all 
functions of r as determined by geometry from the location and pointing direction of the 
telescope). N{D;a,e,i) is the number of particles having effective diameter D and orbits 
with semimajor axis, a, eccentricity, e, and inclination, i. 

We evaluate the integral in Eq. (j6]) by numerical renormalization (see NVBS06). F{D, r) 
is calculated as described in §3.3. N{D] a, e, i) is obtained from numerical simulations in §3.2. 
S{R, L, B) uses theoretical expressions for the spatial distribution of particles with fixed a, 
e and i, and randomized orbital longitudes (Kessler, 1981; NVBS06). 

We assume that the telescope is located at {xt = rtcos0t,l/i = Ttsin0t,-2t = 0) in the 
Sun-centered reference frame with = 1 AU. Its viewing direction is defined by a unit vector 
with components (x^, Uv, ^v)- In Eq- ©, the pointing vector can be also conveniently defined 
by longitude / and latitude b of the pointing direction, where = cos 6 cos/, = cos 6 sin/, 
and Zv = smb. As described in §2, we fix the solar elongation /© = 90° and calculate the 
thermal flux of various particle populations as a function of b and wavelength. The model 
brightness profiles at 12, 25 and 60 /xm are then compared with the mean IRAS profiles 
shown in Fig. El 

To check our code, known as Synthetic InfraRed Telescope (SIRT), we programmed a 
particle version of the algorithm, which should be in many ways similar to the core algorithm 
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in SIMUL (Dermott et al. 1988; see also Dermott et al., 2001). The particle version inputs 
the orbital elements of particles obtained in the orbital simulations (§3.2) and produces their 
orbit clones that are uniformly spread over 27r in mean anomaly M. Thus, every test particle 
is assumed to trace a cloud of real particles having the same orbit as the test particle but 
different angular locations along the orbit. See Vokrouhlicky et al. (2008; their §2.5) for 
a technical description of the algorithm. This procedure is based on the assumption that 
any concentration of particles with a specific M value would be quickly dispersed by the 
Keplerian shear. We employ this procedure to improve the resolution. Without it, the 
number of integrated test particles would be too small to obtain a useful result. 

To be able to compare the results of the particle algorithm with SIRT, the particle 
algorithm must also use smooth distributions in perihelion longitude w and nodal longitude 
Vt. This is achieved by generating additional orbital clones with xu and VL uniformly spread 
over 2tt. Figure H] shows examples of the results obtained from the particle algorithm and 
SIRT. The agreement between the two codes is excellent which gives us confidence that both 
codes work properly. We find that the SIRT algorithm based on the Kessler distribution is 
much faster than the particle one. For this reason, we use the original SIRT code in this 
study. 

4. Results 

Our primary model parameters are the relative contribution of different sources to the 
zodiacal cloud. The total model flux as a function of the ecliptic latitude is obtained as 

Fraodc\{b) = CtAST-^AST + ajFC-^JFC + ttHTC-^HTC + ttOCC-^OCC 5 (7) 

where a are coefficients that satisfy oast + «jfc + c^htc + «occ = 1; and -Fast, -Fjfc, -Frtc 
and -Fqcc are model fluxes obtained for different sources. We normalize them so that the 
ecliptic model flux from each source is equal to that of the mean observed flux at 6 = 0. 
Coefficients a therefore give the relative contribution of different sources at the ecliptic. 

The model flux profiles depend on the particle size, wavelength, and for JFCs also on the 
assumed value of tjFc- As described in §3.2, we tracked particles with 10 < < 1000 /xm. 
These different sizes are treated individually in Eq. ([7]). In particular, we do not attempt to 
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construct plausible size-frequency distributions for different sources. It is therefore assumed 
that a single characteristic particle size, or a narrow range of sizes, can be effectively used 
to model the observed MIR flux. This assumption needs to be verifled later. 

In 4.1, we flrst consider a model where the lifespan of particles is limited by their 
dynamical lifetime. Effects of particle disruptions are discussed in 4.2. 

4.1. Collision- Free Model 

Figure [5] shows the 25-fim flux of -D = 100 fim particles produced by different source 
populations. Note that these proflles do not sensitively depend on the particle size (see Figs. 
E] and [7] for the results for different D). Instead, the main differences between results for 
different source populations in Fig. O reflect the initial orbit distribution of particles in each 
source and their orbit evolution. Therefore, these proflles can help us to identify the source 
population(s) that can best explain IRAS observations. 

The asteroidal particles produce a proflle with a very sharp peak centered at the ecliptic. 
The emission from asteroidal particles near the ecliptic poles is relatively faint. The polar 
emission comes from the particles that evolved by PR drag from a > 2 AU to -R = 1 AU. 
While most asteroidal particles indeed reach 1 AU, they pass too briefly near b = ±90° to 
produce important polar fluxes. This is why most radiation is received from 6 ~ 0, where 
the telescope collects the thermal emission of particles over a wide distance range. A broader 
distribution of orbital inclinations is apparently needed to match IRAS measurements. 

The proflle produced by HTC particles is much broader than the observed one (Fig. 
In this case, the magnitude of the polar fluxes is ^1/2 of that near the ecliptic. This result 
reflects the very broad inclination distribution of HTCs (§3.1). A potentially signiflcant 
contribution of HTCs to the zodiacal cloud is also problematic because the two large HTCs, 
109P/Swift-Tuttle and IP/Halley, tend to librate about mean-motion resonances, causing 
relatively stable orbits for long periods of time. Thus, dust released by HTCs is expected to 
be concentrated along certain location on the sky making it difficult to explain the smooth 
proflle of the zodiacal dust. Note also that Altobelli et al. (2007) have not detected HTC 
particle impacts in the Cassini dust experiment, indicating that HTC particles are relatively 
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sparse. 



The OCC particles, which have a nearly isotropic inclination distribution, produce the 
MIR flux that is constant in latitude (not shown in Fig. [5]). Therefore, the ecliptic and 
polar fluxes from OCC particles are roughly the same and do not match observations. We 
conclude that a single-source model with either asteroidal, HTC or OCC particles cannot 
match the observed profile of the zodiacal cloud. 

We are left with JFCs. It is notable that the width and shape of the JFC profile in 
Fig. E] closely matches observations. The only slight difference is apparent for large ecliptic 
latitudes where the model flux, shown here for D = 100 /im and tjpc = 12, 000 yr, is slightly 
weaker than the one measured by IRAS. We will discuss this small difference below and show 
that it could be explained if: (1) slightly smaller JFC particles were used, and/or (2) the 
zodiacal cloud has a faint isotropic component. We thus believe that the close resemblance 
of our model JFC profile with IRAS data is a strong indication that JFCs are the dominant 
source of particles in the zodiacal cloud. 

Since asteroids and active JFCs have similar inclination distributions (Hahn et al., 2002), 
it may seem surprising that JFC particles produce substantially wider MIR flux profiles 
than asteroidal particles. By analyzing the results of our numerical integrations we find that 
the encounters with terrestrial planets and secular resonances are apparently not powerful 
enough to significantly affect the inclination distribution of drifting asteroidal particles. The 
inclination distributions of the asteroidal particles and their source main-belt asteroids are 
therefore essentially the same (^10° mean i). On the other hand, we find that JFC particles 
are scattered by Jupiter before they are able to orbitally decouple from the planet and drift 
down to 1 AU. This results in a situation where the inclination distribution of JFC particles 
is significantly broader (^20° mean i for R < 3 AU) than that of their source JFCs. This 
explains Fig. [5] and shows limitations of the arguments about the zodiacal cloud origin based 
on the comparative analysis of sources (e.g., Hahn et al., 2002). 

We will now address the question of how the MIR fluxes from the JFC particles depend 
on D and tjFc- We define: 
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where Fjras is the mean IRAS flux, a is the standard deviation of Fjras determined from 
the spread of representative IRAS scans for each h (§2), and -Fmodei = -^jfc (i-e., ajFC = 1 
and Oast = ttHTC = Q^occ = in Eq. ([7])). Note that the integration in Eq. ([8]) is set to 
avoid the intervals in h with strong galactic emission. 

While the deflnition of rf in Eq. ([H]) is similar to the usual statistic (e.g., Press et 
al., 1992), we will not assign a rigorous probabilistic meaning to the rf values obtained from 
Eq. (IH]). This is mainly because it is not clear whether the a values computed in §2 from 
the IRAS data can adequately represent the measurement errors. Instead, we will use Eq. 
([8]) only as an indication of whether a particular model is more reasonable than another one. 
Models with rf ^ 1 will be given priority. For a reference, the JFC model in Fig. [5] gives 
rf = 5.1. 

We calculated rj"^ as a function of D and tjpc to determine which values of these pa- 
rameters flt IRAS observations best. We found that the best flts with rj"^ < 10 occur for 
30<D< 100 fim and tjpc < 30, 000 yr. 

Figure [7] illustrates how the shape of the 25-;um profile produced by JFC particles 
depends on D and tjpc- The profiles become wider with increasing D and tjpc values. For 
^JFC = 12, 000 yr, the best results were obtained with D = 30 fim and D = 100 fim {rf = 3.2 
and 5.1, respectively). The profiles for D = 10 fim are too narrow and clearly do not fit data 
well [r]'^ = 70), while those for D = 1000 /im are slightly too wide (?7^ = 58). We also found 
that there are no really good fits with tjpc > 30, 000 yr, because the profiles are too broad 
near the ecliptic independently of the particle size. 

The best single-source fits discussed above have rj"^ > 1 which is not ideal. According 
to our additional tests this is probably not due to the coarse resolution and studied range of 
D and tjpc- Instead, this may point to: (1) a subtle problem with our JFC model, and/or 
(2) the possibility that additional minor sources (such as asteroids, OCCs or HTCs) should 
be included in the model. Option (1) is difficult to test unless a better model of the JFC 
population becomes available. Here we concentrate on (2) because an ample evidence exists 
(e.g., for ~10% near-ecliptic asteroid contribution from asteroid dust band modeling) that 
these additional sources may be relatively important. 

We start by discussing the results obtained by assuming that the zodiacal cloud has two 
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sources. The motivation for considering the two-source model was the following. First, we 
wanted see whether a combination of two sources could successfully fit the observed profile. 
Second, we attempted to place upper limits on the contributions of asteroid, OCC and HTC 
sources. While it is obvious that models with more than two sources can be tuned to fit 
the data better, it is not clear whether more than two sources are actually required. Our 
two-source models were used to test these issues. 

In the first test, we used the two-source model with asteroids and OCCs (i.e., ajpc = 
ttHTC = in Eq. ([7])). We found that this particular model produces unsatisfactory fits 
ijf > 10) to IRAS observations for all particle sizes considered here (10-1000 /im) (Fig. [8^). 
The model profile is significantly narrower near the ecliptic, where the asteroid component 
prevails, and is too wide for |6| > 40°, where the OCC component prevails. This happens 
mainly due to the fact that the asteroid dust is confined to the ecliptic plane and produces 
a very narrow profile near the ecliptic (Fig. [5l). We also find it unlikely that two so distinct 
populations of objects, such as the main belt asteroids and OCCs, would have comparable 
dust production rates. Thus, we believe that the two-source model with asteroid and OCC 
dust can be dismissed. 

In the second test, we set aocc = «htc = and considered models of the zodiacal 
cloud with the JFC and asteroid components. We found that a small contribution of asteroid 
dust can improve the fits. For example, rf = 0.92 for the D = SO-fim JFC particles with 
«JFC = 0.9 and tjpc = 12,000 yr, and D = 100-/im asteroidal particles with oast = 0.1. 
This represents a significant improvement from t]"^ = 3.2 that we obtained for a single- 
source model with JFC particles only. Values oast ^ 0.3 are clearly unsatisfactory because 
1]"^ > 10 for a AST > 0.3. Also, rj'^ > 3.1 for oast > 0.2 which shows that the fit does not 
improve when we add a >20% asteroid contribution. These results suggests that a very large 
asteroid contribution to the zodiacal cloud can be ruled out. This agrees with the conclusions 
of NVBS06 who found that q;ast ^ 0.15 from modeling of the main asteroid dust bands. 

Finally, the two-component model with the JFC and isotropic OCC sources can fit data 
very well (Fig. [8]d). With D = 100 fim and tjpc = 12,000 yr, corresponding to one of our 
best single-source fits with JFCs, ajFc = 0.97 and aocc = 0.03, we find that r]"^ = 0.36, by 
far the best fit obtained with any two-source model. As Fig. [Hb shows, the fit is excellent. 
We may thus find an evidence for a small contribution of OCC particles to the zodiacal 
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cloud. A much larger OCC contribution is not supported by the data because the fit gets 
significantly worse for aocc > 0.1. For example, rf > 10 for aocc > 0.15 which is clearly 
unsatisfactory. A large contribution of OCC particles can therefore be rejected. 

We propose based on the results described above that the zodiacal cloud has a large JFC 
component (ajpc ^ 0.9), and small asteroid/OCC components (oast ^ 0.2 and aocc ^ 
0.1). To verify this conclusion, we considered three-component models with q;htc = 
and used ajpc, oast and aocc as free parameters. We found that the best two-source 
fit with ajFc = 0.97 and aocc = 0.03 cannot be significantly improved by including a 
small asteroid contribution. Similarly, the fit with ajpc = 0.9 and oast = 0.1 cannot 
be improved by adding a small OCC contribution. Thus, the asteroid/OCC contributions 
cannot be constrained independently because their effects on the combined profiles can be 
compensated by adjusting the D and tjpc values of the dominant JFC particles. 

If we set the parameters of the dominant JFC particles to be D = 100 /im and tjFc = 
12,000 yr, however, the oast and aocc values can be constrained much better (Fig. [9]). 
For example, models with rf < 1 require that oast < 0.22 and aocc < 0.13. Thus, under 
reasonable assumptions, the contribution of asteroid particles to the near-ecliptic IRAS fiuxes 
is probably < 10-20%, in agreement with the results obtained in NVBS06 from modeling of 
the asteroid dust bands. This means that asteroid dust contributes only by <10% to the 
overall zodiacal dust emission at MIR wavelengths. The thermal emission of OCC particles 
can constitute as much as ~10% of the near-ecliptic emission with k.'o% providing the best 
fits (Fig. [9]). When integrated over latitude, the overall OCC component in the zodiacal 
cloud is likely to be <10%. 

For the sake of consistency with the results suggested from modeling of the asteroid 
dust bands (see §3.1; NVBS06), we impose a small asteroid contribution in the JFC/OCC 
model. Figure [lO] shows our preferred fit at different IRAS wavelengths. The rf values of 
this fit in different wavelengths are the following: 0.29 at 12 /im, 0.35 at 25 /im and 0.06 
at 60 yum. This is very satisfactory. Since our model does not include detailed emissivity 
properties of dust grains at different wavelengths (§3.3), we set the emissivity at 25 /im to 
be 1 and fit for the emissivities at 12 and 60 /im. We found that the relative emissivities at 
12 and 60 /im that match the data best are 0.76 and 0.87, respectively. Such a variability of 
MIR emissivity values at different wavelengths is expected for D 100 /im silicate particles 
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with some carbon content (NVBS06). Note also that our preferred D values {D ^ 100 //m) 
are within the range of dominant sizes of particles at 1 AU as determined from spacecraft 
impact experiments {D — 20-200 /xm; Griin et al., 1985). 

4.2. Effect of Disruptive Collisions 

The observational evidence for collisional disruption of interplanetary particles is un- 
deniable (see, e.g., Griin et al., 1985), yet it is very difficult to model the full collisional 
cascade in a computer code as each disrupted dust grain produces numerous fragments. 
The exponentially increasing number of particle fragments, which in reality exceeds 10^^ for 
D > 30 iim (NVBS06), renders any full A'"-body integration impossible. To circumvent this 
problem, the A'"-body integration of a smaller number of "tracer" particles can be coupled 
with a Monte-Carlo model for coUisions as in NVBS06. This method is not ideal. Also, 
any model for the collisional cascade would suffer from our lack of detailed understanding of 
particle properties and their fragmentation during impacts. 

Here wc opt for a very simple approximation of the effect of disruptive collisions. We 
assume that the collisional lifetime of particles is tco\{D) and stop the A^-body integration of 
diameter D particles when t = tcoi{D). Thus, particles keep the same D for t < tcoi(-D) and 
vanish at t = tcoi(-D). This is very a crude approximation of the real collisional cascade, in 
which particles can be eroded by small collisions and do not vanish upon disruptive impacts 
(but produce a range of new particle sizes). Also, tcoi(-D) should be a function of R while 
we assume here that it is not. Still, as we show below, our simple model should be able to 
capture the main effects of particle collisions. 

Our choice of tcoi(-D) is motivated by the published estimates of the collisional lifetime 
of particles based on satellite impact rates and meteor observations. For example, Griin et 
al. (1985) argued that the collisional lifetime oi D = 1 mm particles at 2.5 AU should be 
~10'^ yr (sec also Jacchia and Whipple, 1961). This relatively short lifetime is a consequence 
of the dominant population of D ^ 100-300 fim particles in the inner solar system (e.g.. 
Love and Brownlee, 1993) that are capable of disrupting mm-sized particles upon impacts. 

For comparison, the approximate PR drag timescale of particles to spiral down from 
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2.5 AU to 1 AU is ^pr = 2500 yr x ps, which for p = 2 g cm~^ and s — 500 //m gives 
— 2.5 X 10^ yr. Thus, the PR drag hfetime of these large particles is significantly 
longer than t^oh indicating that they must disrupt before they can significantly evolve by 
PR drag. Using this assumption in the model we found that the profiles produced by large 
JFC particles with tcoi ~ 10^ yr are much narrower in latitude than the ones we obtained 
in §4.1. This is because large particles die before they can evolve to ~ 1 AU, where 
they could contribute to polar fiuxes. The zodiacal cloud cross-section therefore cannot be 
dominated by large JFC grains. The large grains are important to explain radar and optical 
observations of meteors (see §5.3; Taylor and Elford, 1998; Jenniskens, 2006, Wiegert et al., 
2009). 

On the other side of the size spectrum, D < /im particles have tp^ <S ^coi due to 
the lack of small D < 1 /im impactor particles that are blown out of the solar system by 
radiation pressure, and because the PR drag timescale is short for small D (see, e.g., Dermott 
et al., 2001). Thus, the small particles are expected to spiral down by PR drag from their 
initial orbits to it! < 1 AU without being disrupted. Our original results described in §4.1 
are therefore correct for small particles. We showed in §4.1 that D < 30 iim JFC particles 
do not fit IRAS observations well. 

Since tpR -C tcoi for small particles and tp-R ^ icoi for large ones, there must exist 
an intermediate particle size for which tpp^ ~ tcoi- These intermediate-size particles are 
expected to be very abundant in the zodiacal cloud simply because they have the longest 
hfetimes. Griin et al. (1985) and others argued that the transition from the PR-drag to 
collision-dominated regimes must happen near 100 //m. This is consistent with the LDEF 
measurements which imply that the D ^ 200- //m particles represent the dominant mass 
fraction at i? = 1 AU (Love and Brownlee, 1993). 

The question is therefore how to model collisional effects for D ~ 100 fim. This is 
not a simple problem because the effects of the full collisional cascade, including gradual 
erosion of particles and their supply from breakups of the large ones, should be particularly 
important in this transition regime. It would be incorrect, for example, to take the Criin et 
al. estimates of icoi at their face value and remove particles when t > tcoi- In reality, each 
particle can accumulate PR drift during previous stages of evolution when it is still attached 
to its (slightly) larger precursor particles. 
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To test these issues, we assumed a wide range of effective tcoi and calculated model JFC 
profiles for each of these cases. Fig. [TT] shows that the profiles obtained with tcoi < 5 x 10^ yr 
are significantly narrower than the IRAS profiles, even if we tried to compensate for part 
of the apparent discrepancy by OCC particles (Fig. [TTb). On the other hand, profiles with 
tco\ > 6 X 10^ yr are almost indistinguishable from the original results that we obtained 
in §4.1 with tcoi = oo. The transition between 5 x 10^ yr and 6 x 10^ yr occurs near the 
mean PR-drag lifetime of = 100 /xm JFC particles in our model. It clearly makes a larger 
difference whether particles are allowed to drift down to i? = 1 AU or not. The main lesson 
we learn from this exercise is that IRAS observations imply that the zodiacal cloud particles 
have been significantly affected by PR drag. 



4.3. Additional Constraints 

Additional constraints on the micrometeoroid environment near 1 AU are provided by 
radar and optical observations of meteors. For example. Hunt et al. (2004) determined 
the meteor entry speeds from the high-gain ALTAIR radar. For 30 km s^^, the minimum 
detectable mass is 10~^ g (corresponding to -D = 100 fim for p = 2 g cm"'^), while only 
mm-sized and larger meteoroids can be detected by the ALTAIR radar for <20 km s~^. The 
ALTAIR measurements represent a significant improvement in sensitivity relative to that of 
previous radar programs (e.g., Taylor, 1995; Taylor and Elford, 1998). For example, Taylor 
(1995) cited the minimum detectable mass of 10~^ g at 30 km s^^ for the Harvard meteor 
radar, corresponding to -D > 500 fim particles. 

In §5.3, we estimate that the mean atmospheric entry speed of D ~ 100 fim zodiacal 
cloud particles is ~14 km s~^, and that >90% impact at <20 km s^^. Thus, the ALTAIR 
and Harvard radars cannot detect the bulk of small zodiacal cloud particles impacting Earth 
at low speeds. These measurements are instead sensitive to large meteoroids, which carry 
relatively little total mass and cross-section, have short tcoi, and are expected to impact on 
JFC-like orbits. This explains why radar observation show little evidence for populations 
of small particles with orbits strongly affected by PR drag. In Fig. [121 we compare the 
atmospheric entry speeds of D = 1-mm JFC particles with tcoi = 10^ yr with the Harvard 
radar data. This figure documents the dominant role of large JFC particles in meteor radar 
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observations. 

The spatial distribution of sporadic meteors shows several concentrations on the sky 
known as the helion/anti-helion, north/south apex and north/south toroidal sources (e.g., 
Younger et al., 2009, and the references therein). Wiegert et al. (2009) have developed a 
dynamical model to explain these observations. Their main result concerns the prominent 
helion/anti-helion sources for which the large JFC particles clearly provide the best match. 
Our model for large JFC particles is in many ways similar to that of Wiegert et al. (2009). It 
should therefore be consistent with the observed relative strength of the helion/anti-helion 
sources. The more recent high-gain antenna observations show that smaller meteoroids 
appear to show a weaker helion/anti-helion source of eccentric short-period orbits (Mathews 
et al., 2001; Hunt et al, 2004; Galligan and Baggaley, 2004). This imphes that orbits of 
smaller particles should be more affected by PR drag (in agreement with §4.2). 

The motion of interplanetary particles can be probed by high-resolution spectral obser- 
vations of the zodiacal cloud. Reynolds et al. (2004) measured the profile of the scattered 
solar Mg I A5184 Fraunhofer line in the zodiacal cloud. The measurements indicate a sig- 
nificant population of dust on eccentric and/or inclined orbits. In particular, the inferred 
inclination distribution is broad extending up to about 30°-40° with respect to the eclip- 
tic plane. The absence of pronounced asymmetries in the shape of the profiles limits the 
retrograde population of particles to less than 10% of the prograde population. 

These results are in a broad agreement with our model. As we discussed in §4.1, small 
JFC particles are scattered by Jupiter before they are able to orbitally decouple from the 
planet and drift down to 1 AU. This results in a situation where the inclination distribution 
of JFC particles is broad and extends beyond 20°. The model eccentricities of JFC particles 
show a broad range of values with most having e = 0.1-0.5 (see §5.3). This is in a nice 
agreement with the analysis of Ipatov et al. (2008) who found that e ~ 0.3 best fits the 
Reynolds et al. data. 
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5. Implications 

Given the results described in §4 we are now in the position to estimate the total cross- 
section and mass of particles in the zodiacal cloud, the current and historical accretion rates 
of dust by planets and the Moon, and discuss the implications of our work for studies of 
micrometeorites and debris disks. We address these issues below. 

5.1. Zodiacal Cloud Mass 

According to our preferred model with the dominant contribution of JFC particles to 
the zodiacal cloud, the inner circumsolar dust complex has the total cross-section area of 
(2.0 ± 0.5) X 10^^ km^. This is a factor of ~10 larger than the cross-section of asteroidal 
particles in the main asteroid dust bands (NVBS06). The uncertainty of the total cross- 
section was estimated from the range of values obtained for models with rf < 3. Also about 
40% of the total estimated cross section of the zodiacal cloud, or faS x lO^'^ km^, is in particles 
that reside inside Jupiter's orbit (i.e., with R < 5 AU). 

The estimated values are comparable to the effective emitting area of the zodiacal cloud 
defined as 1 ZODY in Gaidos (1999) (1 ZODY = 1.7 x 10^^ km^ assuming blackbody emission 
at 260 K and a bolometric luminosity of 8 x 10~^ L0, where Lq is the Sun's value; Reach 
et al., 1996; Good et al., 1986). Note, however, that we estimate in §5.5 that the bolometric 
luminosity of the inner zodiacal cloud is ~2.5 times larger than the one assumed by Gaidos 
(1999). 

The total mass of the zodiacal cloud is a function of the unknown particle density 
and loosely constrained dominant particle size. With p — 2 g cm~^ and D — 200 //m, we 
estimate that the total mass is mzoDY = 5.2 x 10^^ g, which is roughly equivalent to that 
of a 37-km-diameter body. The zodiacal cloud therefore currently contains relatively little 
mass. Note that these estimates apply to the inner part of the circumsolar dust complex 
that is detectable at MIR wavelengths. The outer circumsolar dust complex beyond Jupiter 
is hkely more massive due to the contribution from KB particles (e.g., Landgraf et al., 2002; 
Moro-Martm and Malhotra, 2003). According to Greaves et al. (2004), the KB dust disk 
may represent up to ~1.2 x 10^^ g. This is up to ~3 x 10^ times the mass of the inner 
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zodiacal cloud estimated here. Note that this is an upper bound only and that the real KB 
dust disk can be much less massive. 

Our mass estimate is at least a factor of ~2 uncertain. For example, if p = 1 g cm~'^ or 
D = 100 fim, we find that mzoDY = 2.6 x 10^^ g. These values are a factor of ~2-4 lower 
than the mass of the zodiacal cloud suggested by NVBS06 from modeling of the asteroid 
dust bands. NVBS06 assumed that the radial distribution of zodiacal particles is similar 
to that of the asteroid dust bands, which is incorrect if JFCs are the dominant source. On 
the other hand, NVBS06 determined the realistic size distribution of zodiacal particles by 
tracking the collisional evolution, while we used the single-size distributions here. 

We estimate that >80% of the total mass at <5 AU should be contained in JFC particles. 
Since these particles can efficiently decouple from Jupiter by PR drag, a large fraction of 
the total mass is distributed relatively close to the Sun. [For reference, we find that 53%, 
19% and 3.7% oi D = 10, 100 and 1000 /im particles released by JFCs, respectively, can 
decouple from Jupiter.] Figure [13] shows the mass fraction of JFC particles contained in a 
sphere of radius R around the Sun. The distribution is steep for i? < 5 AU and shallower 
for R > 5 AU refiecting the orbital properties of our model JFC population. About 30% of 
JFC particles, or about 1.6 x 10^^ g in total mass (for p = 2 g cm~^ and D = 200 pm), are 
located within R<A AU. Also, ^10%, or about 5.2 x 10^^ g, has < 2 AU. 

For a comparison, assuming that the asteroidal particles with D = 200 pm and p = 
2 g cm~'^ contribute by 15% to the near-ecliptic MIR fiuxes, we find that the total masses 
in asteroidal particles with R < 2 AU and R < A AU are 5.3 x lO^'' g and 1.3 x 10^^ g, 
respectively. Thus, the total mass (or number) ratio of JFC to asteroidal particles in the 
inner solar system is <10. Note that this estimate applies as far as the size distributions 
of JFC and asteroidal particles in the zodiacal cloud are similar, which is expected because 
both particle population live in the common collisional environment and have similar PR 
drag timescales. 
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5.2. Mciss Influx on Earth 

We used the Opik algorithm (Opik, 1951; Wetherill, 1967) to estimate the terrestrial 
accretion rate of JFC particles expected from our model. For 30 < D < 300 //m, the average 
impact probabihty of JFC particles on the Earth is ~5 x 10~^ yr~^ per one particle in 
the zodiacal cloud. A similar value is obtained if the impact probability is estimated from 
the number of direct impacts recorded by the A'"-body integrator. Thus, in a steady state 
with ~2 X 10^^ g in the zodiacal cloud, we estimate that ~10^ tons of JFC particles should 
be accreted by the Earth annually. This is larger than the nominal Earth's accretion rate 
of 20,000-60,000 tons yr-^ as determined from LDEF (Love and Brownlee, 1993) and the 
antarctic micrometeorite record (Taylor et al., 1996). 

This may imply that the real Earth's accretion rate is somewhat larger than the LDEF 
values. Alternatively, the LDEF constraints may imply that the real mass of the zodiacal 
cloud is lower than the one estimated here. As we discussed above, the mass of the zodiacal 
cloud estimated here from the IRAS data is at least a factor of ~2 uncertain. It is thus 
plausible that mzoDY ~ 10^^ g (e.g., if p = 1 g cm~^), which would bring our results into 
a better agreement with LDEF. Additional uncertainty in these estimates is related to the 
effects of coUisional disruption of particles and continuous size distribution. 

For comparison, if we assumed that D = 200- /im asteroidal particles are producing 
the full near-ecliptic MIR flux measured by IRAS, the estimated terrestrial accretion rate 
of asteroidal particles would be 1.5 x 10^ tons yr~^. According to NVBS06 and the results 
obtained here, however, the asteroidal particles contribute by only ~10% of the near-ecliptic 
MIR flux. Thus, we find that the asteroid particle accretion rate should be ~15,000 tons yr~^, 
or only ~15% of the JFC particle accretion rate. The asteroidal particles should therefore 
represent a relatively minor fraction of IDPs and micrometeorites in our collections. This 
explains paucity of the ordinary chondritic material in the analyzed samples (see, e.g., Genge, 
2006). 

Using the same assumptions, we estimate from our model that 16,000 tons yr~^ and 
1,600 tons yr^-*^ of JFC particles should be accreted by Mars and the Moon, respectively. 
The accretion rate of JFC particles on the Moon is thus only about ~2% of the Earth's 
accretion rate. This corresponds to a smaller physical cross-section and smaller focusing 
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factor of the Moon. The mass influx on Mars is ~20% of the Earth's accretion rate. For a 
comparison, 1,600 tons and 100 tons of asteroidal particles are expected to fall on Mars and 
the Moon annually. 

Love and Brownlee (1993) found from the LDEF impact record that D ^ 200 fim 
particles should carry most of the mass of zodiacal particles near 1 AU, while we find here 
that D 100 /im provides the best fit to IRAS observations. This slight difference may be 
related to some of the limitations of our model. It can also be real, however, because the 
LDEF size distribution computed by Love and Brownlee (1993) is bending from the steep 
slope a,t D > 300 //m to shallow slope at < 50 //m. The cross-section area should therefore 
be dominated by smaller particles than the mass. From Fig. 4 in Love and Brownlee (1993) 
we estimate that D 100 //m particles should indeed dominate the total cross-section area 
of the zodiacal cloud at 1 AU. 

5.3. Micrometeorites 

These results have implications for the origin of micrometeorites (MMs). MMs are 
usually classified according to the extent of atmospheric heating they endure (e.g., Engrand 
and Maurette, 1998). Cosmic spherules are fully melted objects. Scoriaceous micrometeorites 
are unmelted but thermally metamorphosed objects. The fine-grained MMs and coarse- 
grained MMs are unmelted objects which can be distinguished on the basis of their grain 
size. Based on bulk composition, carbon content, and the composition of isolated olivine 
and pyroxene grains, fine-grained micrometeorites and scoriaceous MMs, which appear to be 
thermally metamorphosed fine-grained micrometeorites, are likely related to carbonaceous 
chondrites. It has been estimated that the ratio of carbonaceous to ordinary chondrite 
MMs is ~6:1 or larger (see, e.g., Levison et al., 2009). This stands in stark contrast to the 
terrestrial meteorite collection, which is dominated by ordinary chondrites. 

A possible solution to this discrepancy is that a large fraction of the collected microm- 
eteorites are particles from the Jupiter-family comets. This possibility has to be seriously 
considered because we find here that the carbonaceous JFC grains should prevail, by a large 
factor, in the terrestrial accretion rate of micrometeoroids. It has been suggested in the past 
that a possible problem with this solution is that the cometary particles should encounter 
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the Earth at large velocities (e.g., Flynn 1995), so that they either burn up in the atmo- 
sphere or are converted into cosmic spherules. Thus, while cometary particles could produce 
fully melted objects such as the cosmic spherules it was not clear whether the less thermally 
processed carbonaceous MMs, such as the fine-grained and scoriaceous MMs, may represent 
cometary material. 

By assuming tcoi ^ 5 x 10^ yr for D = 100 /im as required by IRAS observations (see 
§4.2), we find from our model that the mean impact speed of D ~ 200 fim JFC particles on 
Earth is ^14.5 km s~^ (Fig- [HI see §4.3 for a discussion of the size dependence of impact 
speed and its relevance to meteor observations). This value is only slightly higher than 
that of the asteroidal particles (^12.5 km s^^). The comparable impact speeds of JFC and 
asteroidal particles in our model are a consequence of PR drag which efficiently circularizes 
the orbits before they can reach 1 AU (Fig. [T5l) . We thus find that the impact speeds of the 
JFC particles are low and do not pose a serious problem. Based on this result and the high 
terrestrial accretion rate of JFC particles on Earth (§5.2), we propose that the carbonaceous 
MMs in our collections are grains from the Jupiter-family comets. A large contribution from 
primitive material that may have been embedded into the main asteroid belt according to 
Levison et al. (2009) is probably not needed. 



5.4. Historical Brightness 

It is believed that the main source of JFCs is the scattered trans-Neptunian disk, which 
should have decayed by a factor of ~100 over the past ~4 Gyr (LD97; Dones et al., 2004). 
If the JFC population decayed proportionally, we can estimate that the ecliptic component 
of the zodiacal dust should have been ~100 times brighter initially that it is now. This 
corresponds to the near-ecliptic 25- /im fiux of about 7 x 10^ MJy sr~^. 

A different insight into the historical brightness of the zodiacal cloud can be obtained in 
the framework of the Nice model (Tsiganis et al., 2005), which is the most complete model 
currently available for the early evolution of the outer solar system. In the Nice model, 
the giant planets are assumed to have formed in a compact configuration (all were located 
at 5-18 AU). Slow migration was induced in these planets by gravitational interaction with 
planetesimals leaking out of a massive primordial trans-planetary disk. After a long period 
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of time, most likely some 700 Myr after formation of the giant planets (Gomes et al., 2005), 
planets crossed a major mean motion resonance. This event triggered a global instability that 
led to a violent reorganization of the outer solar system. Uranus and Neptune penetrated 
the trans-planetary disk, scattering its inhabitants throughout the solar system. Finally, the 
interaction between the ice giants and the planetesimals damped the orbits of these planets, 
leading them to evolve onto nearly circular orbits at their current locations. 

The Nice model is compelling because it can explain many of the characteristics of the 
outer solar system, (Tsiganis et al., 2005; MorbideUi et al., 2005; Nesvorny et al., 2007; 
Levison et al., 2008; Nesvorny and Vokrouhhcky, 2009). In addition, the Nice model can 
also provide an explanation for the Late Heavy Bombardment (LHB) of the Moon (Tera 
et al., 1974; Chapman et al., 2007) because the scattered inhabitants of the planetesimal 
disk, and main belt asteroids destabilized by planetary migration, would provide prodigious 
numbers of impactors in the inner solar system (Levison et al., 2001; Gomes et al., 2005). 

Assuming that the historical brightness of the zodiacal cloud was proportional to the 
number of primitive objects that were scattered into the inner solar system on JFC-like 
orbits, we can estimate how it changed over time. In the pre-LHB stage in the Nice model, 
the leakage rate from the planetesimal disk beyond 15 AU was likely not significant relative 
to that at LHB. We thus expect that the MIR emission from the inner zodiacal cloud at 
R < 5 AU should have been relatively faint, except if a massive population of particles 
was sustained by collisions in the pre-LHB asteroid belt. Here we focus on the LHB and 
post-LHB stages. 

According to Wyatt et al. (2007), the astcroidal debris disk is expected to decay by 
orders of magnitude from the time of Jupiter's formation, which marked the start of the 
fragmentation-dominated regime in the asteroid belt (e.g., Bottkc et al., 2005), to LHB. It 
thus seems unlikely that a massive population of debris could be sustained over 700 Myr 
by the collisional grinding of main belt asteroids. Instead, it has been suggested that the 
collisional grinding in the massive trans-planetary disk at i? > 15 AU should have produced 
strong MIR emission peaking at ~100 /xm (Booth et al., 2009; hereafter B09). Being more 
distant the trans-planetary disk probably decayed more slowly by collisions than the asteroid 
belt. Thus, in the pre-LHB stage, the Wien side of the trans-planetary disk emission may 
have exceeded the one from the inner zodiacal cloud down to ~20 //m (B09). 
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During the LHB, as defined by the Nice model, large numbers of outer disk planetes- 
imals were scattered into the inner solar system and the inner zodiacal cloud could have 
become orders of magnitude brighter than it is now. To estimate how bright it actually was, 
we used simulations of the Nice model from Nesvorny and Vokrouhlicky (2009; hereafter 
NV09). NV09 numerically tracked the orbital evolution of 4 outer planets and 27029 objects 
representing the outer planetesimal disk. The mass of the disk was set to be 35 Earth masses. 
In total, NV09 performed 90 different numerical integrations of the Nice model, only some 
of which ended with the correct orbits of the outer planets. 

We used these successful simulations to determine the number of scattered objects with 
JFC-like orbits as a fraction of the total initial number of planetesimals in the trans-planetary 
disk. Figure [16] shows how this fraction changed over time in one of the NV09 successful 
simulations. Consistently with the estimated physical lifetime of modern JFC (LD97) we 
assumed that the physical lifetime of planetesimals after reaching g < 2.5 AU for the first 
time was 10^ years. Objects past their physical lifetime did not contribute to the statistic. 

Immediately after the planetary instability occurred in the Nice model, the estimated 
fraction of planetesimals having JFC-like orbits was ~7 x 10~^ (Fig. fTB]) . It then decayed 
to ^10~^ at 50 Myr after the start of LHB. Even though the NV09 simulations gradually 
loose resolution at later times due to the insufficient number of tracked particles, we can 
still estimate that the fraction was ~10~^ at 500 Myr, or about 3.4 Gyr ago in absolute 
chronology. 

Charnoz et al. (2009) and Morbidelli et al. (2009) argued, using the crater record 
on lapetus and the current size distribution of Jupiter's Trojans, that the total number of 
D > 2 km planetesimals in the pre-LHB trans-planetary disk was ~10^'^-10^^. Using this 
value and Fig. [161 we find that there were ~7 x 10^ JFCs with D > 2 km at time of the 
LHB peak, tLHB, and ~10^ JFCs at tLHB + 50 Myr. These estimates are at least an order of 
magnitude uncertain mainly due to the poorly known size distribution of small planetesimals 
in the trans-planetary disk. 

For comparison, Di Sisto et al. (2009) found, in a good agreement with the previous 
estimates of LD97, that there are ^100 JFCs with D > 2 km and q < 2.5 AU in the current 
solar system (with about a factor of 50% uncertainty in this value). Therefore, if the inner 
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zodiacal cloud brightness reflects variations in the size of the historical JFC population, we 
find that it has been ~7 x 10'^ brighter at tLHB and 2 x 10^ brighter at tLHB + 50 Myr 
than it is now. This would correspond to the near-ecliptic 25-/im fluxes of 5 x 10® and 
10^ MJy sr"^, respectively. These values largely exceed those expected from dust particles 
that were scattered from the trans-planetary disk (B09). Most of the action was apparently 
over by tLHB + 500 Myr, when our model suggests that the inner zodiacal cloud was only 
~10 times brighter than it is nowE 



5.5. Distant Observations of the Zodiacal Cloud 

Figure [I7| shows how the present zodiacal cloud would look like for a distant observer. 
If seen from the side, the brightest inner part of the zodiacal cloud has a disk-like shape with 
a ~1.6 ratio between the ecliptic and polar dimensions. Similar shapes have been reported 
by Hahn et al. (2002) from Clementine observations of scattered light. At a larger distance 
from the Sun, the shape of the zodiacal cloud resembles that of a walnut. The axial ratio 
becomes ^1.3 at i? = 5 AU. 

The radial brightness profiles in Fig. [17] show a steep dimming of the zodiacal cloud 
with R. For R < 1 AU, a factor of ~10 in brightness is lost per 1 AU. For 1 < i? < 5 AU, 
factor ~10 is lost per 2 AU. These profiles are approximate because we ignored the effect 
of collisions in our model, which should be especially important for -R < 1 AU. It is unclear 
how the shape of the zodiacal cloud would look for _R > 5 AU because we did not model the 
contribution from KB dust. 

Figure [18] shows the Spectral Energy Distribution (SED) for distant unresolved obser- 
vations of the zodiacal cloud. At a distance of 10 pc from the Sun, SED of the present inner 
zodiacal cloud is 1.4 x 10~^ Jy at 24 fim and 5.5 x 10"^ Jy at 70 fim, corresponding to 
the excesses over the Sun's photospheric emission at these wavelengths of about 3.4 x 10~^ 
and 1.1 X 10"'^, respectively. For comparison, the approximate 3a excess detection limits 



^ These estimates should only be taken as a rough guideline to the historical zodiacal cloud brightness 
because the collisional environment in the dense disk of JFC particles at LHB must have been very different 
from the one existing today. It is therefore not exactly correct to assume that the historical brightness of 
the zodiacal cloud was strictly proportional to the population of JFCs. 
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of Spitzer telescope observations of Sun-like stars are 0.054 at 24 fim and 0.55 at 70 /xm 
(Carpenter et al., 2009). The MIR emission of the present inner zodiacal cloud is therefore 
undetectable by distant unresolved observations with a Spitzer-class telescope. Specifically, 
the detectable emission levels are ^160 and ~500 larger at 24 and 70 /xm, respectively, than 
those of the present inner zodiacal cloud. 

When the flux is integrated over wavelengths, we find that the fractional bolometric 
luminosity of the inner zodiacal cloud, Lzody, relative to that of the Sun, Lq = 3.839 x 
10^6 W, is Lzoby/Lq ~ 2 X 10"''. This is a larger value than 10"®-10~'' suggested by 
Dermott et al. (2002) and perhaps comparable to that of KB dust at > 30 AU (Stern et 
al. 1996). The effective blackbody temperature of the zodiacal cloud can be estimated from 
Teff = 5100/Amax, wherc A^ax is the wavelength of the SED maximum in microns. With 
Amax = 18.5 fim, this gives this gives Tes = 276 K. 

5.6. MIR Excess during LHB 

B09 studied how the MIR excess of the solar system debris disk varied with time. 
According to them, the main source of the pre-LHB MIR emission should have been the 
population of dust particles produced by collisions in the massive trans-planetary disk at 
i? > 20 AU. In Fig. [181 we show the model SED produced by the B09 trans-planetary disk. 
Being dominated by collisions (as opposed to PR-drag regime; see Wyatt, 2005), the trans- 
planetary particles are destroyed before they could evolve to i? < 20 AU. The SED emission 
therefore peaks at longer wavelengths (~100 /im) than the SED of the present zodiacal cloud 
(~20 /im). Also, with ~35 Earth masses in the pre-LHB trans-planetary disk, its estimated 
MIR emission is strong and produces excesses of ~0.1 at 24 fim and ~50 at 70 /im over the 
Sun's photospheric emission at these wavelengths. These values are comparable to those of 
observed exozodiacal debris disks (B09). 

The trans-planetary disk objects, including small dust particles, became scattered all 
around the solar system during LHB. This led to a significant depletion of the trans-planetary 
particle population which could not have been compensated by the collisional cascade be- 
cause collisions became increasingly rare in the depleted disk. The MIR excess should have 
thus dropped by orders of magnitude within several hundred Myr after the LHB start. B09 
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estimated that the 24-/im excess of dust particles scattered from the trans-planetary disk 
should have dropped to ~3 x 10^^ at the present time. This is about an order of magnitude 
lower value than the 24 /im excess estimated by us for the current inner zodiacal cloud. 
Thus, there must have been a transition epoch some time after LHB when the 24 /im ex- 
cess stopped being dominated by dust particles scattered from the trans-planetary disk and 
became dominated by particles produced by JFCs. 

Figure [19] illustrates how the MIR emission of the inner zodiacal cloud should have 
varied with time during LHB. The size of the JFC population was estimated by using 
the methods described in §5.4. We then scaled up the MIR emission of the present in- 
ner zodiacal cloud by the appropriate factor (see §5.4). We found that the 24 fim excess 
reached FzoDY/-fsun(24/im) ~ 20 at the LHB peak and stayed for about 100 Myr above 
the Spitzer's 3a detection limit. It dropped down to ~10 times the value of the present 
zodiacal cloud at tLHB + 500 Myr. We were unable to determine how this trend continues 
after tLHB + 500 Myr because of the resolution issues with the NV09 simulations (§5.4). We 
expect that FzoDY/-fsun(24/im) should have decayed by an additional factor of ~10 from 
to + 500 Myr to the present time. 

The 70- /xm excess behaves similarly (Fig. [T9b). It reaches FzoDY/-^Sun(70/im) ~ 70 
at the LHB peak and decays at later times. Given the tighter detection limit of Spitzer at 
70 fim, the 70-yum excess would remain detectable by Spitzer for 50 Myr, which is roughly 
half of the interval during which the 24-/im excess could be detected. Also, when these values 
are compared to the ones estimated by B09 for the pre-LHB trans-planetary disk, we find 
that the 70-/im excess expected for JFC particles at the LHB peak is comparable to that 
of the pre-LHB excess produced by collisions in the trans-planetary disk. We would thus 
expect that the solar system's LHB has not produced any significant increase of the 70-/im 
emission. 

Conversely, the 24-/xm excess raises a factor of ~100 above the B09 pre-LHB level, 
indicating that the solar system became significantly brighter during LHB at these shorter 
wavelengths. During LHB, the emission from JFC particles should have exceeded that of the 
scattered trans-planetary particles by ~20 (compare Fig. [T9b with Fig. 5 in B09). Thus, 
the system does not need to have a significant cold dust disk at the same time as the hot 
dust disk to provide material for the hot disk. In fact, we find here that the hot disk can 
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be fed by D >1 km objects, which have httle total cross-section area to be detected in the 
cold disk, but large mass to sustain the hot disk upon their disintegration at <10 AU (see 
§6). Since the decay rates of both populations after LHB should have been similar, their 
ratio should have remained roughly constant over time suggesting that the trans-planetary 
dust did not represent any significant contribution to the post-LHB emission of the zodiacal 
cloud at 24 iim. 

5.7. Debris Disks 

These results have interesting implications for our understanding of hot debris disks ob- 
served within 10 AU around mature stars (TriUing et al., 2008; Meyer et al., 2008; Carpenter 
et al., 2009; see also a review by Wyatt 2008). It has been argued that some the observed 
brightest hot disks, such as HD 69830, HD 72905, HD 23514, rj Corvi and BD+20307, can- 
not be explained by assuming that they are produced by the coUisional grinding of the 
local population of asteroids (Wyatt et al., 2007). Specifically, Wyatt et al. (2007) pointed 
out that the emission from a locally produced population of debris is expected to be much 
weaker than the observed emission because disks become depleted over time by collisions. 
This problem cannot be resolved by assuming a more massive initial population because the 
massive population would decay faster. Instead, Wyatt et al. proposed that the bright hot 
debris disks can be seen around stars with planetary systems that are undergoing the LHB 
instability akin to that invoked in the Nice model. 

In §5.6, we estimated how the 24- //m and 70-//m excesses varied during the solar system's 
LHB. We found that the 24-//m excess should have rapidly risen by a large factor from the 
pre-LHB value and then gradually decayed. It would remain detectable by a Spitzer-class 
telescope for about 100 Myr after the LHB start, or ~2% of the Sun's current age. The 
24-//m excess reached values >10 at LHB, which is comparable to those of the brightest 
known hot disks (Wyatt et al., 2007). Conversely, the solar system's LHB has not produced 
a significant increase of the 70-//m excess relative to the pre-LHB level. The 70-//m excess 
decayed after LHB and became undetectable by a Spitzer-class telescope after ~50 Myr. 
Thus, if the timing is right, debris disks may show the 24-/xm excess but not the lO-jim 
excess. This could explain systems such as HD 69830, which shows a large excess emission 
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at 8-35 /xm (Beichmann et al. 2005), but lacks 70-//m emission, and HD 101259 (Trilling et 
al., 2008). 

In a broader context, our study of the zodiacal cloud implies that: (1) the populations of 
small debris particles can be generated by processes that do not involve disruptive collisions 
(see §6); and (2) observed hot dust around mature stars may not be produced from a 
population of objects that is native to <10 AU. Instead, in the solar system, most particles 
located within the orbit of Jupiter are fragments of planetesimals that formed at >15 AU. 
These icy objects are transported to <5 AU by gravitational encounters with the outer 
planets and disintegrate into small particles by disruptive splitting events (thought to occur 
due to processes such as the pressure build-up from heated volatiles or nucleus spin-up; see 
§6). If these processes are common around stars harboring planets, the coUisional paradigm 
in which debris disks are explained by collisions may be not as universal as thought before 
(see, e.g., Wyatt 2008). 

5.8. LHB Accretion Rates 

If our basic assumptions are correct, large quantities of dust should have been accreted 
by the Moon, Earth and other terrestrial planets during LHB. For example, assuming 10^^ g 
y^-i jj^gg^j^ accretion rate over 100 Myr we estimate that ~10^^ g of extraterrestrial material 
should have fallen on the Earth at the time of LHB (with ~50% of this mass accumulating 
in the first 10 Myr). This is ~50 times more mass than the quantity accumulated by the 
Earth at its current accretion rate over 4 Gyr. The Moon should have accreted about 2% of 
the Earth's value during LHB, or ~ 2x10^° g in total over 100 Myr. These estimates are at 
least one order of magnitude uncertain. 

For comparison, the mass of large impactors estimated from the number and size dis- 
tribution of lunar basins is 6 x 10^^ g (Hartmann et al., 2000). Thus, the total mass of 
dust deposited on the Moon during LHB should have been only ~l/30 of that of the large 
impactors. 

The LHB is of fundamental interest in studies of the origin of life because it immediately 
precedes the oldest evidence for a biosphere (Awramik et al., 1983; Schidlowski, 1988; Mojzsis 
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et al., 1996). The significance of our results in this context is that JFC dust grains can bring 
in unaltered primitive material from the outer solar system. They could potentially be the 
source of the earliest organic material that gave rise to hfe on Earth (e.g., Jenniskens, 2001; 
Jenniskens et al., 2004). [Asteroids are an important source of IDPs but they can accrete 
material from as far as ~4 AU. It is not likely that organic material at such distances can 
survive the T Tauri wind of the young Sun.] 

6. Comet Disruptions 

Today, ^3.8 Gyr after LHB, the steady flux of JFCs from the outer solar system is 
keeping the zodiacal cloud at roughly constant brightness. We find from our numerical 
simulations that the mean dynamical lifetime oi D — 200 //m JFC particles is 10^ yr. Thus, 
to keep the zodiacal cloud at constant brightness, a continuous input of ~3.4 x 10^^/10^ = 
3.4 X 10^^ g yr""^, or roughly 1,100 kg s~^ is required in our model. This estimate is robust 
because it is insensitive to the assumed p and D values of particles (i.e., lighter particles have 
shorter dynamical lifetimes). It neglects, however, the loss of particles due to the disruptive 
collisions. The real input rate should therefore be slightly larger, probably somewhere in 
the 1,000-1,500 kg s"^ range. This is only shghtly larger than 600-1000 kg s~^ suggested by 
Leinert et al. (1983) from modeling of the Hehos 1 and 2 data. 

For comparison. Reach et al. (2007) suggested from the Spitzer survey of cometary 
debris trails that the total meteoroid input from active short-period comets is ~300 kg s~^. 
This is ~3-5 times lower value than what would be required, according to our estimate, to 
keep the zodiacal cloud brightness at constant brightness. While some of the uncertainties in 
our model and the Reach et al. results may be blamed for this discrepancy, we believe that 
this comparison may indicate that the trails of active comets represent only a fraction of the 
real mass loss in comets. In fact, it has been suggested that the main mass-loss mechanism 
in comets is their spontaneous (i.e., non-tidal) disruptions followed up by the progressive 
splitting of comet components into smaller fragments (e.g., Weissman, 1980; also see Chen 
and Jewitt, 1994; Boehnhardt, 2004; Fernandez, 2005; Jenniskens, 2006). 

The best documented case of comet fragmentation is that of sungrazers. These are small 
comet fragments that are detected because they pass very close to the Sun and are seen in 
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backscattered light by solar telescopes (Sekanina and Chodas, 2004, 2005). Specific cases of 
JFCs that were observed to spontaneously split or break up into two or more components 
include 5 IP/ Harrington, 73P/Schwassmann-Wachmann 3 and 141P/Machholz 2 (Fernandez, 
2005). Observations of these events show that there does not seem to be a correlation between 
the splitting event and orbital phase of the parent object, which provides motivation for how 
particles were released from JFCs in our model (§3.1). 

Several fragmentation mechanisms may to explain the splitting of cometary nuclei: (1) 
rotational splitting when the centrifugal force exceeds nucleus' self-gravity and material 
strength; (2) sphtting by thermal stress produced by the variable distance to the Sun; and (3) 
splitting by internal gas pressure caused by sublimation of subsurface pockets of volatile ices 
(e.g., CO). It has not been possible find the main culprit so far. Plausibly, several different 
mechanisms contribute and more observational constraints will be needed to distinguish 
between them. See Weissman (1980) and Boehnhardt (2004) for a discussion. 

Fernandez (2005) compiled a list of 12 observed split JFCs. He found that the chance of 
JFC undergoing an observed splitting event is ~1% per orbital period. This should be taken 
as a lower limit on the actual number of splitting events because many are undetected. For 
example, Chen and Jewitt (1994) estimated that a comet has a ~1% chance to split per yr. 
Thus, over its active lifespan of about 10"^ yr (LD97), typical JFC would undergo as many 
as ~100 sphtting events. These events may lead to the situation where the comet nucleus 
becomes completely dissolved into small particles. The zodiacal cloud may thus plausibly be 
sustained by disintegrating Jupiter-family comets. 

Our order-of-magnitude estimate supports this possibility because JFCs evolving into 
the inner solar system represent a continuous input of mass that is apparently large enough 
to compensate for the zodiacal cloud mass loss. Moreover, we found no evidence in this work 
for tjFc values larger than the physical lifetime of active comets estimated in LD97. Most 
JFC comets should therefore be dissolved on timescales comparable to their active lifetime. 

Using the size distribution of JFCs from Tancredi et al. (2006), we find that the total 
mass of JFCs with radius 0.1 < s < 10 km and q < 2.5 AU is 3.9 x 10^^ g. Assuming that 
this mass is injected into the zodiacal cloud every 10^ yr (LD97), we find the total mass input 
of 12,000 kg s~^. This is significantly larger than the mass input required to maintain the 
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zodiacal cloud in a steady state (1,000-1,500 kg s~"^), possibly suggesting a ~10% yield of the 
disintegration process. Note, for example, that some JFCs or their large fragments can be 
removed (e.g., impact planets or leave the solar system) before they could fully disintegrate. 
Also, icy particles released by comets subhmate at < 5 AU and do not contribute to the 
inner zodiacal cloud. 

Di Sisto et al. (2009) determined the physical lifetime of JFCs to be ~3 times shorter 
than LD97. Using Di Sisto et al. estimate, we find that the JFC population should require 
the mass input of 35,000 kg s~^. The yield of the disintegration process may thus be as 
low as ~3%. For comparison, Di Sisto et al. found the following fractions of JFCs that are 
completely dissolved by splitting events: 51% for radius s = 1 km, 13% for s = 5 km and 
8% for s = 10 km. 

The initial size distribution of particles resulting from the splitting process is uncertain, 
but meteor showers from freshly ejected dust trails, such as Phoenicids, indicate that the 
distribution should be fairly flat with most mass in mm to cm size grains. This initial size 
distribution is modifled by collisions as JFC particles decouple from Jupiter and drift to lower 
R where collisions are more common. As discussed in §4.2, the collisional effects explain why 
D fti 100 fim provide the best fit to the IRAS data, because these intermediate-size particles 
have longest lifespans (e.g., Griin et al., 1985; Dermott et al., 2001). 

Additional evidence that disruptions/splitting events of JFCs may dominate the popula- 
tion of interplanetary particles in short-period orbits comes from observations and modeling 
of the meteor showers. Specifically, it has been estabhshed that most meteor streams were 

produced by recent (<few thousand years ago) comet disruptions (see Jenniskens 2008 for 
a review). For example, 1956 Phoenicids and near-Earth object 2003 WY25 are most likely 
fragments produced by a breakup of D 1819 Wl (Blanpain) (Jenniskens and Lyytinen, 2005; 
Watanabe et al., 2006). In addition to 2P/Encke, there are other known comet fragments 
moving in the Taurid stream (Jenniskens, 2006), also pointing to a disruption event. Gem- 
inids. Phaeton, and 2005 UD can also be linked to the common parent body (Jenniskens, 
2006; Ohtsuka, 2005; Jewitt and Hsieh, 2006). The type of disintegration that produced 
these large fragments and meteoroid streams is probably like that of the 1995 breakup of 
73P/Schwassmann-Wachmann 3, which will cause a shower of tau-Herculidis in 2022. 
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The meteoroid streams that were associated with comet disruptions are much stronger 
than the meteoroid streams produced by active JFCs. Thus, the strong meteoroid streams 
may represent an important hnk between JFCs and the zodiacal cloud. They should become 
increasingly more dispersed due to effects of planetary perturbations. Eventually, the parti- 
cles should be well mixed in orbital space, producing both the sporadic meteoroid complex 
and zodiacal cloud. Notably, the time- integrated flux of visual meteors at Earth is dominated 
by about a factor of ~10 by sporadics (Jones and Brown, 1993). 

Based on modeling of meteor radar observations, Wiegert et al. (2009) demonstrated 
that the prominent helion/anti-helion pair of sporadic meteors is most likely produced by 
JFCs. This result provides further support to the zodiacal cloud model proposed in this work 
because it shows that the JFC particles are an important part of the zodiacal dust complex 
at 1 AU. According to Wiegert et al., the north/south apex pair is probably produced by 
retrograde long-period comets, perhaps suggesting an OCC component in the zodiacal cloud. 
As we showed in §4, a small contribution of isotropic OCC particles is also required to explain 
IRAS observations of the zodiacal cloud. 

7. Compcirison with Previous Work 

The origin and evolution of the zodiacal cloud has been the subject of numerous studies. 
For example, Liou et al. (1995) suggested, based on modeling in many ways similar to our 
own, that the observed shape of the zodiacal cloud can be accounted for by a combination 
of ~l/4-l/3 of asteroid dust and ~2/3-3/4 cometary dust. We found a much larger JFC 
contribution and much smaller asteroid contribution in this work. The cause of this difference 
is unknown. Possibly, it stems from some of the approximations used by Liou et al. (1995). 
For example, they used particles from comet 2P/Encke to represent the whole population of 
particles released by JFCs. This comet has a special orbit (a = 2.22 AU and q — 0.33 AU) 
that is not representative for the JFC population as a whole. 

Different constraints on the origin of the zodiacal cloud have been obtained from model- 
ing the asteroid dust bands. For example, Dermott et al. (1994b) suggested that the particles 
originating in the main asteroid belt supply ~l/3 of the zodiacal cloud, while NVBS06 es- 
timated the contribution of asteroidal particles to be <10%. Our results presented in §4 
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are more in line with the NVBS06 estimate. Specifically, we found that a >20% asteroid 
contribution to the near- ecliptic MIR fluxes can be ruled out from IRAS observations. If 
correct, this limits the asteroid contribution to the overall cross-section of the zodiacal cloud 
to a sub-10% level. 

Hahn et al. (2002) used Clementine observations of the zodiacal cloud at optical wave- 
lengths and arguments based on the inclination distribution of small bodies in the solar sys- 
tem to argue that at least ^90% of the zodiacal cloud cross section enclosed by a 1-AU-radius 
sphere around the Sun is of cometary origin. They also found that ^45% optical cross-section 
at 1 AU comes from JFCs and/or asteroids. Unfortunately, a distinction between JFC and 
asteroid dust could not have been made because Hahn et al. used an approximate model for 
the interplanetary dust complex. According to our model, the contribution of JFC is much 
larger than the one found by Hahn et al. (2002). Thus, while we agree with the general 
conclusion of Hahn et al. about the predominant comet dust population, our results are 
more specifically pointing out JFCs as the main source. 

8. Origin of Particle Populations beyond Jupiter 

Our findings are in a broad agreement with the results obtained from dust detectors 
onboard spacecrafts. For example, Altobelli et al. (2007) identified two main groups of 
particles in the Cassini's Cosmic Dust Analyzer data set (measurements in the ecliptic plane 
between Jupiter and Saturn). The first group of impactors consists of particles on bound 
and prograde orbits, most probably having moderately eccentric and moderately inclined 
orbits. These grains are consistent with JFCs. Impactors of the second group were identified 
as small interstellar dust particles, perhaps including a minority of beta meteoroids. 

Landgraf et al. (2002) reported results from the dust experiments onboard the Pioneer 
10 and 11 spacecrafts. They found that the spatial number density of >10 fim particles at the 
ecliptic is only slowly declining with heliocentric distance in the 3-18 AU range. Specifically, 
there is no obvious gap beyond 4 AU, expected if asteroidal particles were the dominant 
source of dust in the inner solar system (Fig. [T3|l . 

The nearly constant spatial density of the circumsolar dust beyond 5 AU is puzzling. To 
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explain it, Landgraf et al. (2002) proposed that particle populations beyond Saturn are be 
dominated by dust produced in KB collisions (see also Moro-Martm and Malhotra, 2003). 
The observed radial density profile beyond 5 AU is produced in their model by combining 
the contributions from KB particles, whose spatial density raises with R, and cometary 
particles, whose density dechnes with R. Indeed, the spatial density of JFC particles that 
we obtain from our model rather steeply declines with R at R > 5 AU. Thus, the Kuiper 
belt dust may indeed be needed to explain Pioneer measurements. [A possible caveat of 
these considerations is that the impact rates measured by Pioneer 1 and 2 should be mainly 
those of ~10 fim particles, while the dominant size of particles in the inner zodiacal cloud 
is ~100-200 /xm.] 

An alternative possibility is that we do not correctly determine the distribution of JFC 
grains for i? > 5 AU in our model. This alternative is attractive for the following reasons. 

If the trans-Neptunian population is in the collisional equilibrium for D < 1 km, most 
mass should be contained in comet-size and larger bodies rather than in D < 1000 fim grains. 
Since the transfer of this material to the Jupiter-crossing orbit is sizc-indcpcndent (driven 
mainly by the encounters to outer planets), JFCs must represent much more mass than their 
grain-sized orbital counterparts. Thus, assuming that JFCs can be efficiently dissolved by 
splitting events, the dust population they produce should be much more important than the 
one evolving from the Kuiper belt in the form of dust grains. 

Di Sisto et al. (2009) found a very high sphtting rate of JFCs with only a shallow 
dependence on their perihehon distance (oc q" with a ~ —0.5). Thus, if JFCs can be 
efficiently dissolved at large q, the radial distribution of JFC dust should significantly differ 
from the one obtained here (see §3.1 for our assumptions). The spatial density of JFCs is 
proportional to R^ with 7 ~ 0.5 (LD97; Di Sisto et al., 2009). Since |a| ~ 7, the number 
of splitting events, and therefore the number of generated JFC particles, should be roughly 
independent of R. It might thus be plausible to explain the Pioneer measurements with JFC 
particles alone, i.e., without a major contribution from KB particles. A detail investigation 
into these issues goes beyond the scope of this paper. 
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9. Summary 

We developed models for various source populations of asteroid and cometary dust 
particles. These models were based on our current understanding of the origin and evolution 
of asteroids, Jupiter-family, Halley-type, and Oort-Cloud comets. We launched sub-mm 
particles from these populations and tracked their orbital evolution due to radiation pressure, 
PR drag and planetary perturbations. The thermal MIR emission from the synthetic particle 
distributions were determined and the results were compared to IRAS observations. 

The main goal of this modeling effort was to determine the relative contribution of 
asteroid and cometary material to the zodiacal cloud. We found that asteroidal particles 
produced by the main belt collisions cannot produce the zodiacal cloud emission at large 
ecliptic latitudes simply because the main belt asteroids have generally small orbital incli- 
nations, and because the orbital effects of planetary encounters and secular resonances at 
a < 2 AU are not powerful enough to spread the asteroid dust to very large orbital inclina- 
tions. Therefore, most MIR emission from particles produced in the asteroid belt is confined 
to within ~30° of the ecliptic (Fig. [6]). Conversely, the zodiacal cloud has a broad latitudinal 
distribution so that strong thermal emission is observed even in the direction to the ecliptic 
poles (Fig. E]). 

Based on the results discussed in §4, we proposed that >90% of the zodiacal cloud 
emission at MIR wavelengths comes from dust grains released by Jupiter-family comets, 
and ^10% comes from the Oort cloud comets and/or asteroid collisions. We argued that 
disruptions/splitting events of JFCs are more likely to produce the bulk of observed dust in 
the inner solar system than the normal JFC activity. The relative importance of JFC and 
Kuiper-belt particles beyond Jupiter has yet to be established. 

Using our model results, we estimated the total cross-section area and mass of particles 
in the zodiacal cloud, current and historical accretion rates of dust by planets and the Moon, 
and discussed the implications of our work for studies of micrometeorites and debris disks. 
We found that JFC particles should dominate the terrestrial accretion rate of micromete- 
oroids. This may explain why most antarctic micrometeorites have primitive carbonaceous 
composition. If the spontaneous comet disruptions are also common in the hot exozodiacal 
debris disks, the coUisional paradigm used to explain their properties may not be as universal 
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as thought before. 
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Fig. 1. — The upper solid lines in each panel show IRAS scan 180_24 (see Table 3 in NVBS06) 
that has been smoothed by a low pass-filter to remove point sources and instrumental noise. 
Different panels show fluxes at 12, 25, 60 and 100 iim IRAS wavelengths. The bottom 
solid lines show the contribution of three main asteroid dust bands to the observed fluxes. 
According to NVBS06, these dust bands contribute to the observed fluxes by f=:i9-15% within 
10° to the ecliptic, and <5% overall. The strong signal at 100 /im between latitudes b ~ —80° 
and b w —30° is the galactic plane emission (also apparent at 60 /im). Figure from NVBS06. 
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Fig. 2. — Mean IRAS profiles at 12, 25 and 60 fim wavelengths. To make these profiles, 
the selected IRAS scans were centered at the ecliptic, smoothed by a low-pass filter, and 
combined together. The gray rectangles at I < —78° and 40° < I < 70° block the latitude 
range where the mean fluxes were signiflcantly affected by the galactic plane emission. We 
do not use the excluded range in this work. The uncertainties of the mean flux values are not 
shown here for clarity; they are too small to clearly appear in the plot. The characteristic 
errors at different wavelengths averaged over latitudes are (7i2|um = 0.59 MJy sr~^, cr25^m = 1-1 
MJy sr~^ and (Teo/im = 2.7 MJy sr^^. They increase with wavelength due to the larger role 
of galactic emission at longer wavelengths. 
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Fig. 3. — Orbital distribution of JFCs from the LD97 model. Panels show the perihelion 
distance (a), and inclination (b), as functions of the semimajor axis for JFCs with q < 2.5 
AU and tjFc = 12,000 yr. See §3.1 for the definition of tjFC- The 2:1 and 3:2 mean motion 
resonances with Jupiter correspond to the gaps in the distribution at a ^ 3.3 and 3.96 AU, 
respectively. The inclination distribution of JFCs shown here is remarkably similar to that 
obtained by Di Sisto et al. (2009; their Fig. 10). 
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Fig. 4. — Comparison between results obtained from the particle algorithm (dots) and SIRT 
(solid lines). Case 1 corresponds to asteroidal particles with a = 2.5 AU, e = 0.1 and i = 10°. 
Case 2 corresponds to cometary particles crossing the Earth's orbit with a—1 AU, e = 0.5 
and i — 50°. In both cases we assumed that particles have D — 100 fim and are distributed 
randomly in Jl, zu and M. The flux at 25 /im was normalized to a population of 10^^ particles 
in Case 1 and 2 x 10^^ particles in Case 2. Observations with — 1 AU and Iq — 90° vere 
assumed. In Case 2, the particle algorithm shows a scatter around the exact solution due to 
the rough resolution of the distribution near the telescope's location. We used 5 x 10^° orbit 
clones in the particle algorithm. 
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Fig. 5. — Comparison of the 25-/im profiles produced by different sources with IRAS obser- 
vations. The black line shows our mean IRAS scan for Iq = 90°. The colored lines show 
profiles expected from different source populations: asteroids (green), JFCs (red) and HTCs 
(blue). The OCC flux, not shown here for clarity, is a nearly constant function of latitude. 
The maximum flux in each profile has been normalized to 1. We used D = 100 fim and 
tjFc = 12, 000 yr. The main differences between profiles are not sensitive to the exact choice 
of D and other model parameters. 
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Fig. 6. — Dependence of the shape of 25-//m profiles produced by asteroidal particles on 
D. The dashed hne shows the mean 25-//m IRAS profile for Iq — 90°. The upper solid 

curves show the model results for the same wavelength and elongation. The bottom lines 
show the residual flux obtained by substracting the model fiux from the mean IRAS profile. 
Results for D = 30, 100 and 300 /xm asteroidal particles are shown with slightly broader 
profiles corresponding to larger D. The profiles for = 10 and 1000 fim, not shown here, 
are narrower than the ones for D = 30 /^.m. For D = 1000 /xm, this is mainly due to the 
effects of disruptive collisions that destroy large grains before they could evolve down to 1 
AU (see discussion in §4.2). None of the model profiles obtained with asteroidal particles 
can match IRAS observations. 
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Fig. 7. — Dependence of the shape of 25-//m profiles produced by JFC particles on D and 
^JFC- The panels show results for different tjpc: (a) ijpc = 12, 000 yr, (b) ^jpc — 30, 000 yr, 
(c) ijFC = 50, 000 yr and (d) tjpc = 100, 000 yr. The dashed hne in each panel shows the 
mean 25-//m IRAS profile for Iq = 90°. The upper solid curves show the model results for 
the same elongation. The bottom lines show the residual fiux obtained by substracting the 
model fiuxes from the mean IRAS profile. Results for D = 10, 30, 100, 300 and 1000 fim are 
shown in each panel with broader profiles corresponding to larger D. Some of these model 
profiles do not match IRAS observations well. Specifically, tjpc > 50, 000 yr, D > 300 //m 
and D = 10 /im can be clearly ruled out. 
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Fig. 8. — Examples of fits where we modeled the zodiacal cloud as having two sources. Fluxes 
at 25 fim are shown, (a) The best-fit model with asteroid and OCC sources. This model 
does not fit IRAS observations well. The model profile is too narrow near the ecliptic and 
too wide overall, (b) Our best two-source model. Here we used q;jfc = 0.97, aocc = 0.03, 
D = 100 iim and tjpc = 12, 000 yr. The faint isotropic component improves the fit quality 
so that rj"^ — 0.36 in (b). This may suggest that the zodiacal cloud contains a small but 
significant fraction of OCCs particles. 
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Fig. 9. — Model constraints on the contribution of asteroid and OCC particles to the zodiacal 
cloud. Here we used a three-source model with D = 100 /im and ohtc = 0. For a range of 
the ttocc and oast values, we set ajFc = 1 — «occ — «ast, and calculated r]'^ (Eq. [H]) for 
each model. The contours show rj'^ = 3, 10 and 30. The shaded area denotes the parameters 
of our best-fit models with r]'^ < 1. These models have aocc < 0.13 and aocc < 0.22 thus 
placing an upper limit on the near-ecliptic contribution of asteroid and OCC particles. 
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Fig. 10. — Our preferred fit with ctjFc = 0.85, aocc = 0.05, q;ast = 0.1 and q;htc — 0. 
Particles with D = 100 //m and ^jpc = 12, 000 yr were used here. The dashed fines sfiow tfie 
mean IRAS profiles at 12, 25 and 60 //m. The upper and lower sofid fines are the model and 
residual profiles, respectively. The wiggle in the residual profiles for |6| < 10° may occur due 
to a slight problem with our asteroid dust band model. 
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Fig. 11. — Dependence of the shape of 25-/im profiles produced hy D = 100 fim JFC particles 
on tcoi- The dashed line in each panel shows the mean 25-yLim IRAS profile for Iq = 90°. 
The upper sohd curves show the model results for the same elongation. The bottom fines 
show the residual flux obtained by substracting model fluxes from the mean IRAS proflle. 
Results for icoi — 10^, 10^, 5 x 10^ and 10^ yr are shown in each panel with broader proflles 
corresponding to the larger tcoi values. In panel (a), we show results for the single-source 
model with JFC particles only. The results of the two-source model with JFC and OCC 
particles are illustrated in (b). We included the OCC component in the model to try to 
compensate for the defficient polar fluxes from JFC particles with short tcoi- Profiles with 
^coi ^ 5 X 10^ yr do not match IRAS observations well. 
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Fig. 12. — Comparison of atmospheric entry speeds oi D — 1-mm JFC particles with icoi = 
10^ yr with the Harvard meteor radar data (Taylor, 1995). There is a good agreement 
between the two distributions for >20 km s~^. The number of impacts from large JFC 
particles drops at <15 km s~^. The Harvard data is affected by strong biases for <20 km 
s~^, because the detectable ionization level produced by a meteor is a strong function of 
meteor speed. 
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Fig. 13. — Cumulative distribution of JFC (solid line) and asteroidal (dashed) particles as a 
function of heliocentric distance R. For each R, the value on the Y axis gives the fraction of 
particles (or equivalently fraction of the total mass) contained within a sphere of radius R 
around the Sun. The JFC particles show a shallower slope with about 70% having R > 4 AU. 
Conversely, 99% of asteroidal particles have R < 4 AU. Note that the distributions shown 
here have been normalized to 1 and do not reflect the actual relative contribution of JFC 
and asteroidal particles to the zodiacal cloud. This figure merely shows the trends in both 
populations with heliocentric distance. 
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Fig. 14. — Model distributions of Earth-impact speed of JFC (solid line) and asteroidal 
(dashed) particles with D = 200 fim. Since the effects of the gravitational focusing have 
been accounted for in the calculation, the minimum impact speed is equal to the escape 
velocity from the Earth's surface, or about 11.2 km s~^. Majority of JFC particles have the 
impact speeds in the 11.2-15 km s~^ range. JFC particles with larger impact speeds have 
lower impact probability but are important for interpretation of the meteor radar data (e.g., 
Wiegert et al, 2009). 
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Fig. 15. — Eccentricity (top panel) and inclination (bottom) distributions of JFC particles 
in our model. The dashed hnes show the distributions for all JFC particles with i? < 10 AU. 
The sohd lines show the distribution for 0.9 < it! < 1.1 AU. The upper plot illustrates that 
the orbits of JFC particles drifting by PR drag become nearly circularized before reaching 
1 AU. The inclination distribution does not change much during this evolution. 
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Fig. 16. — Number of objects on JFC-like orbits during LHB as a fraction of the total 
number of planetesimals in the pre-LHB trans-planetary planetesimal disk. The fraction 
was determined from the n22 simulation of the Nice model in NV09. We extracted all orbits 
from that simulation with perihelion diatance g < 2.5 AU, orbital period P < 20 yr and 
assumed that the physical lifetime of these objects was 10^ yr (LD97). We also used an 
averaging window of 1 My to improve the statistics. The total mass of the JFC population 
can be estimated from this plot by multiplying the fraction shown here by the initial mass 
of the trans-planetary disk. With the 35 Earth-mass disk, the peak in the mass of the JFC 
population at i ~ corresponds to ~0.3 lunar masses. 



-73- 






Fig. 17. — Zodiacal cloud brightness at 24 //m as seen by an observer at 10 pc. Two 
projections are shown: (top) polar view for an observer with Z = 10 pc; and (bottom) 
side view of an observer in the ecliptic plane (Y = 10 pc). The three isophotes in each of 
the two left panels correspond to 5 x 10"'^, 5 x 10"^ and 5 x 10~^ Jy AU"^ with 1 AU^ 
at 10 pc corresponding to 0.01 arcsec^. The shading scale is linear in log^g of brightness. 
The right panels show the brightness variation with the heliocentric distance along the cuts 
denoted by the dashed lines in the left panels. There are two lines in the bottom-right panel 
corresponding to the polar and ecliptic profiles. 
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Fig. 18. — Spectral density distribution of the present inner zodiacal cloud as seen by an 
observer at distance 10 pc form the Sun. For reference, we also plot SED of the Sun and the 
pre-LHB trans-planetary disk as determined by Booth et al. (2009). The two arrows show 
the approximate 3(7 detection limits of the Spitzer telescope at 24 and 70 /xm (Carpenter et 
al., 2009; Wyatt et al, 2008). 
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Fig. 19. — Expected variation of excesses at 24 /xm (panel a) and 70 /xm (b) during LHB 
(solid lines). To determine the excess values at different times during LHB, we used Fig. [16] 
to estimate the number of objects that were scattered from the trans-planetary disk into 
the JFC-like orbits. By comparing this number to the present population of JFCs, a scale 
factor has been determined to represent the brightness increase of the inner zodiacal cloud 
over its current value. The discontinuity in the lines near t = 30 Myr appears because we 
changed the size of the averaging running window, 5t. For t < 30 Myr, we used 5t = 1 
Myr; for t > 30 Myr, we used 5t = 50 Myr. The large 6t value is needed for t > 30 Myr to 
improve the statistics. For reference, the plot also shows the values predicted by Booth et 
al. (2009) for the pre-LHB trans-planetary disk, approximate Spitzer detection limits and 
present inner zodiacal cloud (dashed lines). 



